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ABSTRACT 


A simulation model for the open ocean submarine versus 
submarine search and detection problem is presented. The 
objective of the simulation is to estimate the probability 
with which a nuclear powered attack submarine will achieve 
sonar detection of a nuclear powered transiting submarine 
using a search Sian based on external intelligence. A 


detailed description of the model and its use are included 


along with a typical analysis. 
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I. INTRODUCTION 


This thesis contains a detailed description of the 
VECTORED INTERCEPT MODEL (VIM), an event sequenced Monte 
Carlo simulation of a submarine versus submarine open ocean 
search and detection problem assisted by external intelli- 
gence. The objective of VIM i225 estimate the probability 
with which a nuclear powered attack submarine will achieve 
a direct path sonar detection of a nuclear powered target 
submarine based on a given geographical and tactical situa- 
tion. The method used is to estimate the probability of 
detection at various times throughout the trial, combine 
the estimates to form a cumulative trial probability, then 
average over all trials. The method of estimation is 
explained in detail in Chapter III, Section A. 

VIM was initially developed to determine the probability 
of an attack submarine establishing direct path sonar 
contact of a target submarine as a result of meveiwaiue 
intelligence provided by a shore based detection system 
for an analysis of the U.S. strategic defense plan by the 
Joint Chiefs of Staff War Gaming Agency (JWGA) [1]. In 
response to an interest expressed by Commander, Submarine 
Development Group II, VIM was extended into a more general 
form to be used in planning of exercises, extension of 


exercise results, and paremetric studies of the open ocean 
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search and detection system. No reports of use of VIM 
by Submarine Development Group II have been published. 

The twofold purpose of this thesis is to provide a 
detailed description of VIM including instructions in use of 
the model and to illustrate its use by means of an example. 
Chapter I is intended to acquaint the reader with the nature 
of the simulation and general format of the situation being 
modelled. 

Chapter II will assist the reader in understanding the 
operation of VIM and determining if it is suitable for use 
in the analysis of his problem. 

Chapter III contains the details of the simulation 
including the alternatives available to the user in relating 
the model to a specific problem. A description of the over- 
all logic flow and each of the subprograms is included. 

Chapter IV describes the input and savour neenantens 
of VIM including the various options designed to simplify 
data input and display, and an explanation of output 
messages. . 

| Chapter V is a user's manual which includes the defini- 
tions of all input parameters and the keypunch instructions 
for constructing the input data deck. The user's manual is 
cross-referenced to the descriptive passages in the thesis 
which pertain to the inputs, and to the definitions of 
related inputs. 

Chapter VI is a discussion of the methodology and data 


requirements for the execution of a simulation experiment 





using VIM, and an example of a suevied to which VIM was 
applied. 

Appendix A describes a method of utilizing layer depth 
information in an experiment. Appendix B is a program 


Listing. 


A. PROGRAM INFORMATION 

VIM is written in FORTRAN IV for the IBM system 360 
computer. The program space requirements are 120,000 bytes 
for storage and execution. No plotting, card punching or 
special tape drives are employed. Execution time varies 


\ ‘ 
between 0.15 seconds and 2.5 seconds per trial depending 


on the period of simulation and the level of output detail. 





_-II. MODEL DESCRIPTION 


This chapter includes a brief description of the simula- 
tion, an introduction to the terms which will be used 
throughout the thesis, and a discussion of the categories 
of input data. It is intended that within this chapter the 
reader will be able to assess the applicability of VIM to 
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his problem. 


A. MODEL STRUCTURE 

In simulating the situation outlined in Chapter Tt, an 
attempt was made to maintain Simplicity while permitting 
sufficient flexibility so that a wide variety of problems 
could be addressed. The following is a brief synopsis in 
terms of ships, playing area, ship motion, intelligence, 
detection and convergence zones of the method by which VIM 
portrays the basic scenario. 

Ships: Two nuclear submarines are considered. An 
attack submarine is assigned the mission of intercepting 
and establishing close sonar contact within weapon range of 
a target submarine. The mission of the target is to complete 
an assigned transit undetected by the attack submarine. It 

- must be noted that at no time during the simulation does 
the attacker actually establish close contact. Rather, 
# the probability of the event is estimated at well defined 


times throughout each trial. 
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Playing area: The geography of the open ocean is 
approximated by a plane surface of unlimited dimension, 
and with a coordinate system centered at the attacker 
initial position. Ship positions are recorded in terms of 
miles east or west of the origin (X coordinate), and miles 
north or south of the origin (Y coordinate). Operating 
depths are not modelled in VIM, however a method of account-— 
ing for layer depths is presented in Appendix A. 

Ship motion: The attacker is initially assigned to 
a waiting station where it will remain until it receives 
intelligence on the target. This intelligence is received 
either through its own sensors or via a communications link 
with a surveillance facility. “Upon receipt of the intelli- 
gence the attacker attempts to establish sonar contact by 
transitting to a point believed to be ahead of the target 
and then executing a search pattern about that point. The 
attacker remains at the search station until new intelli- 
gence is received. 

Basic target motion is an input to the simulation. The 
target will deviate from its track in a random manner if 
desired, and will also execute an evasion pattern if a 
counterdetection occurs, returning to the track upon comple- 
tion of the evasion maneuver. 

Intellignece: Information concerning the target upon 
which the attacker bases its search pattern may be generated 
by a surveillance facility and transmitted via a communica- 


tions link, or derived directly from a convergence zone 
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detection. Intelligence data includes the contact time 
and estimates of course speed and position. 

Detection: VIM treats detection in one of three ways. 
The target responds to a detection of the attacker by 
executing an evasion pattern; the attacker may respond to 
a convergence zone detection by attempting to intercept ' 
the target; or at potential points of detection the prob- 
ability with which the attacker detects the target is 
estimated with neither unit reacting to the event. 

A counterdetection occurs whenever an evasion plan 
has been provided to the target and the signal excess at 
‘the target as determined from the sonar equation [2] is 
positive. Counterdetections are computed deterministically 
based on ship noise, sonar gain and range between units. 
At various times determined by the trial geometry, the | 
probability that the attacker detects the target is esti- 
mated. This probability is computed using sonar equation 
for signal excess [2] under the assumption that the amount 
of reduction in noise level as sound propagates through the 
water is a normally distributed random variable. These 
estimates are combined to provide an overall estimate for 
each trial. See Chapter III, Section B for the details of 
this process. 

Sonar detection opportunities arising for the attacker 
within a convergence zone (CZ) may result in one of two 
actions depending on the nature of the problem. If a C2 


detection would satisfy the requirements of the attacker's 
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mission, the detection probability estimate is combined 


with other probability estimates to form the overall esti- 
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mate for each trial. If a CZ detection would fail to satis- 
fy the attacker's mission, the detection probability 
estimate becomes the probability of a CZ intelligence 
detection at that point. If a CZ intelligence detection 
oceurs the attacker attempts to intercept the target based 


on the CZ intelligence information. 


, 


B. DATA 

The data required for operating VIM falls loosely 
within six categories: environment, identification, 
“options, situation, target track and intelligence. While 
a certain amount of overlap exists among these categories, 
they are treated as distinct data units by the model. 

Environmental Data: This data group contains the 
radiated and self noise curves of both units along with 
the propagation loss curve for the geographical area being 
considered. 


Identification data: This data group contains informa- 





tion to assist the user in distinguishing the output from 
among various experimental runs and identifying the signif- 
icant features of each run. : 

- Option data: This data group contains information of 
an administrative nature which defines the way in which VIM 


» will treat the various aspects of the simulation. Typical 


of the option data is the number of trials to be run, the 





yy 


level of output detail desired, and the method of generating 
intelligence. 

Situation data: This is a general category which 
includes all data not required elsewhere. Typical situation 
entries are attacker operating speeds, initial target posi- 
tion, target evasion parameters and intelligence error 
parameters. 


Target track data: This data defines the target's 


“transit path. Included with each leg of the target track 


is an estimate of the vulnerability of the target to detec- 

tion by the surveillance facility during the leg. 
Intelligence data: If the observed output of a surveil- 

lance system is to be duplicated, it is entered as intelli- 


gence data. Acquisition time, position error, course and 


speed estimates may be listed. Information not supplied 


will be generated randomly by VIM. 
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IIIT. DETAILED MODEL DESCRIPTION 


This chepter contains a detailed exposition of the 
methods used in VIM to simulate the various facets of the 
submarine versus submarine search and detection problem. 
Section A is a discussion of each phase of the simulation 
including the options available for use in examining a wide 
variety of situations. Section B is an outline of overall 
program flow in terms of trial events. Section C deals 
with the function and logic of each subroutine as it relates 


to the overall program. 


A. MODEL DETAILS 

For the purposes of this discussion, VIM is divided 
into seven categories: target motion, convergence zone 
detections, intelligence, communications, attacker motion, 
detection and evasion. 

Target motion: The target is provided an initial posi- 
tion and a sequence of future courses, speeds and times. 
The target procedes on the first course and speed until the 
problem time indicates a change is required. Except for the 
initial position the target track is independent of target 
location, being defined in terms of course, speed and time 
to start subsequent legs. The target will deviate from its 
assigned track in order to execute an evasion maneuver, but 
will resume the base course and speed when evasion has been 


terminated. If desired, target initial position and the 
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course and speed on each leg can be varied randomly. In 
this case, the course and speed taken by the target after 
evasion will be randomly determined at that time. 

Convergence zone detections: One of two alternative 
interpretations may be selected for CZ detection opportu- 
nities. If a CZ detection of the target by the attacker 
may be adequate to fulfill the attacker's mission, then 
for each period the attacker spends within the convergence 
zone, a Single detection’probability is estimated when the 
range between the ships is as near the central CZ radius as 
it is expected to get. This probability estimate is treated 
in the same way as all other such estimates. If a CZ detec- 
tion would fail to satisfy the attacker's mission the 
detection probability estimate generated within the con- 
vergence zone is understood to be the probability of a CZ 
intelligence detection to which the attacker will respond 
as if it had received intelligence from the surveillance 
facility. CZ intelligence data is generated randomly based 
on convergence zone parameters and cannot be combined with 
surveillance intelligence. CZ detection probability esti- 
mates are not combined with other detection probabilities 
if the CZ intelligence option is exercised. 

Intelligence: Intelligence data generated by a surveil- 
lance facility may be developed in one of several ways: 


all data input; some data input, the remiander created 


within the program; or, all data created within the program. 








Within each of the latter two categories several options 
are available. If all data is to be input, the acquisition 
time of the intelligence, the position error relative to 
the actual target coordinates and the course and speed 
estimates for each surveillence facility detection are 
required. The course and speed estimates may be omitted 
from the input data and will be selected randomly or com- 
puted from position estimates. Figure 1 shows an example 
of typical course and speed estimates based on intelligence 
positions. Omission of the position error from the input 
data will result in a random error generated at the time of 
the intelligence acquisition. If all intelligence is to be 
generated within the program, no intelligence data is 
entered. The times of acquisition are random based on the 
susceptibility of the target to the surveillance system on. 
each leg. Position, course and speed errors are determined 
as indicated above. 

Communications: At periodic intervais throughout each 
trial the attacker monitors an intelligence broadcast which 
provides the information generated by the most recent 
surveillance detection subject to a fixed delay time. The 
interval between broadcasts is an input parameter for the 
attacker transit phase, but is set to two hours during the 
search phase. A continuous broadcast monitoring option is 


available whereby intelligence is provided to the attacker 


immediately following the fixed delay. 








PN 


Figure 1. Intelligence Course and Speed Estimates Based 


on Position Estimates 
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Attacker motion: Attacker motion is governed by one 
of three situations: initial waiting station, response to 
intelligence and search for the target. At the beginning 
of each trial the target is on waiting station proceeding 
on an input course and speed. It will eontinue on the 
initial vector until receipt of an intelligence estimate 
regarding the target. In response to the intelligence, 
transit speed is taken and course is determined to inter- 
cept the target based on’intelligence information. The 
search phase begins when the attacker reaches the expected 
intercept point. The attacker assumes its assigned search 
speed and searches back and forth perpendicular to the 
estimated target tracks on legs of specified duration. The 
search phase continues until new intelligence information 
is received. | 


Detection: The sonar equation for figure of merit [2] 


FM = Nr - max { Ns, Na } + Nrd ~ Nai 


where Nr is the noise level in decibels (db) radiated by 


the source; 


Ns.is the background noise in db at the receiver; 

Na is the ambient noise level of the ocean in db; 

Nrd is recognition differential: the ability of the 
system to distinguish between signal and random 
noise; 

Ndi is directivity index: the ability of the system to 


identify the direction from which a signal is coming; 
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Nrd-Ndi: measured in db is treated herein as the single 


entity, sonar gain. 


Propagation loss, L, is a measure of the reduction in 
noise level in db as sound is propagated through the water. 
L is a function of the range from the noise source and is 
assumed in VIM to have a normally distributed error. 

Signal excess, Se, is defined as Se = FM-L, and is also a 
normally distributed random variable with the same variance 
as L. When Se is zero the probability of detection is 


assumed to be 1/2. 


Detection probability, p, is estimated from the formula 


p = $(Se/s) 


where $(z) is the probability that a stanaard vernal random 
variable will be less than z, and‘s is the standard devia- 
tion of propagation loss. For example, if FM = 109 db, 

L = 95.5 db and s = 9 db, then Se = 13.5 db, Se/s = 1.5 

and p = .93. In developing the aggregate probability esti- 
mate, P, for the trial it is assumed that each estimate is 


independent and 


P = P' + p(1-P') 

where P' was the overall probability prior to the estimate p. 
Probability estimates are generated at the start of each 

event and after any course or speed change by either ship, 


but are only included in the overall estimate under the 


following circumstances: 








1. the attacker is in transit to a search station and 
a course or speed change by either unit results in 


a change from decreasing range to increasing range; 


2. the attacker is in transit to a search station and 


achieves closest point of approach; 


'3. the attacker is on search station and achieves its 
maximum detection probability during that search 


phase; or 


4, the attacker is within a convergence zone and is 
as close as expected to the central radius and the 


CZ intelligence option was not selected. 


Under the CZ intelligence option detection probability 
generated within a convergence zone is taken to be ne. Giobs 
ability of a CZ intelligence detection, and CZ probability 
estimates are not included with the trial estimate. 

Evasion: Four prepared evasion patterns are available 
to the target with provisions for a fifth pattern if desired. 
In pattern 1 the target reverses course and reduces speed 
by one half. Pattern 2 calls fon 3 clockwise course change 
and a opene change. In pattern 3, illustrated in Figure 2, 
the target changes course relative to the bearing of the 
attacker away from attacker track and changes speed. (In 
pattern 4, illustrated in Figure 3, the course change is 
relative to the attacker course away from its track and 
speed is changed. Except in pattern 1 the magnitudes of 


the course and speed changes along with the evasion interval 


ed 





Figure 2. Evasion Pattern 3 
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If relative target bearing is less than 180° the evasion 
angle is subtracted from the reciprocal of target course. 
If relative target bearing is greater than 180° the evasion 


angle is added. 
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Figure 3. Evasion Pattern 4 
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If relative target bearing is less than 180° the evasion 
angle is subtracted from the true bearing of the attacker. 
If relative target bearing is greater than 180° the evasion 


angle is added. 
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are inputs. Evasion is terminated after the fixed interval 
or when a new target track leg is ordered. 

Provision is made for the user to introduce his own 
evasion routine. The calling sequence is controlled by the 
normal input parameters. 

If the target has been assigned an evasion routine, 
counterdetections will occur when a deterministic evaluation 
of the sonar equation indicates Se is positive. Counter- 


detections are supressed during the evasion maneuver. 


B. TRIAL EVENTS 

VIM is composed of a sequence of events, each event 
occuring at a time determined by previous events, randomly 
or input. At each event, the positions of the target and 
attacker are updated, the appropriate action dictated by 
the event is. executed, if necessary new event times are 
generated, data is stored or-printed and the trial clock 
is advanced to the time of the next event. 

Target Course Change: The target assumes bite course 
and speed of its next leg. Upon executing a course change 
the target terminates evasion and sets its course and 
speed to the next assigned values either exactly or randomly. 
Using the new course and speed the expected time of closest 
point of approach (CPA), counterdetection time, detection 
probability and CZ detection time are calculated. 

Counterdetection: A Counterdetection event signals 
either the commencement or termination of evasion. If no 


evasion is specified counterdetections are supressed. At 
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the time of counterdetection the target commences its 
evasion routine and the next counterdetection event is set 
to occur after the evasion interval. Time of CPA, detection 
probability, CZ detection time and, at termination of eva- 
sion, counterdetection time are computed. 

Intelligence Detection: This event generates the 
results of a detection of the target by a surveillance 
facility in accordance with the desired option. The data 
is stored for later release to the attacker. Neither ship 
responds to an intelligence detection. 

Convergence Zone: This event results in a detection 
opportunity for the attacker. If the CZ intelligence op- 
tion has been selected, the detection probability estimate 
is taken to be the probability of a CZ intelligence detec-— 
tion. If a CZ intelligence detection occurs the intelli- 
gence data is computed based on parameters compatible with 
the submarine versus submarine tracking problem, and the 
attacker responds to it as it would to surveillance intelli- 
gence received during a communications period. If a CZ 
intelligence detection fails to occur no action results and 
the event is terminated. When CZ intelligence is not de- 
sired, but a convergence zone is present, the C4 probability 
estimate is combined with other estimates to form the over- 
all trial probability estimate. No response by either unit 


is required. The CZ event occurs once for each passage 


through the convergence zone. 





Communications: This event provides the opportunity 

for new intelligence data from the surveillance facility 
to be provided to the attacker. At the time of a Communi- 
cations event the list of intelligence collection times is 
searched for the latest time which precedes current time by 
more than the fixed delay. If the data has been previously 
transmitted to the attacker or has been superceded by CZ 
intelligence the event is terminated. If the intelligence | 
data has not been previously transmitted a course is deter- 
mined by which the attacker will intercept the estimated tar- 
get. The new CPA time, counterdetection time, convergence 
zone time and detection probability are computed. If the 
attacker is on search station the maximum detection prob- 
ability of that search phase is combined with other proba- 
bility estimates to form the cumulative trial estimate of | 
detection probability. | 

~ Detection Probability: The current detection probabil- 
ity estimate is either combined with other estimates to 
form the overall trial probability or used in the selection 
of a maximum search probability depending on the tactical 
‘situation. “If the attacker is in transit to an intercept 
point, the current probability estimate is combined with 
the previous probability estimates to produce the cumulative 
trial estimate. if the attacker is on search station the 
estimate is compared with all previous estimates generated 
_ during the search pattern and if it is the maximum it is 
. retained for later reference. Otherwise it is dropped and 


the event is terminated. 
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Attacker Course Change: This event indicates an attack- 
er course and speed change upon arrival at search station 
and at the end of each search leg. Course changes related 
to new intelligence are either communications events or 
convergence zone events and are not referred to as attacker 
course change events. The attacker's search station is 
centered at a point determined by the intersection of the 
attacker's transit track and the estimated target track. 
Upon arrival at the center of the search station the attacker 
assumes a course 90 degrees clockwise from the estimated 
target course and procedes at search speed for one half 
the time allotted for a search leg. At the end of the 
search leg the attacker reverses course and searches in 
the new direction for the entire prescribed interval. Each 
attacker course change event is accompanied by computation 
of new CPA time, counterdetection time, detection probabil- 
ity and convergence zone time. 

End of Trial: The trial is terminated when no further 
positive detection probabilities are expected and in any 
case when trial time exceeds a nominal termination time 


by 100 hours. The specific criteria for ending a trial are: 


1. The final intelligence estimate of: the trial has 
been transmitted to the attacker, 

2. The attacker is on search station, 

3. Probability of detection is zero, and 

4, The range of the target to the search station is 


increasing. 
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In addition to the above criteria, if the final intelli- 
gence of the trial results in an estimated speed advantage 
for the target so that intercept appears impossible the 
trial is terminated. 

If the attacker is on search station at the end of the 
trial maximum detection probability generated during the 
final search phase is included in the overall probability 


estimate for the trial. All trial statistics are recorded 


for summary display and the next trial is begun. 


C. LOGIC AND SUBROUTINES 

This section contains a description of the main program, 
MAIN, logic flow, illustrated in Figure 5, along with an | 
explanation of each of the subroutines. Program activities 
can be categorized as run preparation, trial preparation, 
event execution and display, trial completion and run com- 
pletion. Figure 4 shows the relationship of each of the 
categories. The details of each activity are best found by 
examining Figure 5. and the explanations of subroutines 
referenced within the figure. 

MAIN controls the sequence of operations during the 
execution of a series of runs. Run preparations involve 
the reading and display of input data and are carried out 
at the beginning of each run. The environmental inputs 
are read in before the first run only. Trial preparations 
precede each trial and include initializing key variables 


and developing the event list. For each trial, MAIN 
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Figure 4. VIM Logic Summary 
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Figure 5. Main Program Flow Chart 
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Figure 5. (Continued) 
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Figure 5. (Continued) 


Event type 1: Target course change 













Cease evasion 
set new target 
course and 
speed. 


D 


To new event times 


7 Event type 2: Surveillance system detection 








To end of 
event 


Continuous 
broadcast? 








Set the time 
for the next 
communications 
event 


32 





Figure 5. (Continued) 
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Figure 5. (Continued) 


Event type 5: Attacker course change 
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Figure 5. (Continued) 


Event type 7: End of trial 
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Figure 5. (Continued) 
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Figure 5. (Continued) 
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executes the events in chronological order until the trial 
termination criteria are satisfied. Following each trial, 
MAIN directs the data assembly for the trial and the desired 
output, then initiates the next trial or run summary as 
appropriate. Following the run summary, the pattern is 
repeated starting with run preparations until all runs have 
peen completed. 

While primary control of VIM rests within MAIN many of 
the details of the simulation are contained in subroutines 
which are called by MAIN and other subroutines. In the 
text that follows, these subroutines will be: referred to by 
name and are listed in order of their appearance in the 
program listing, Appendix B. 

INOUT: The input and output functions are performed 
primarily by INOUT. These functions are described in 
detail in Chapter IV. In addition to reading all input 
data and writing most of the output messages INOUT initial- 
izes trial variables and prepares the statistical data for 
display. Eight entry points are provided and are called 
only by MAIN at various critical times throughout the seriés 
of puns. . 

Entry INOUT: This routine reads tille eard and 
environmental data is the preliminary step in a series of 
runs. 

Entry INOUT2 : This routine reads run designator, 


run description, options, situation, target track and 
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intelligence data in preparation for each run. Optimal 
tareet iInvercept is computed. 

Entry INOUT3: This section of INOUT displays the 
input agra at the start of each run. 

Entry INOUT4: New random intelligence times are 
printed by this routine at the start of each trial if 
desired. Counters used in compiling run statistics are 
set to zero. 

Entry INOUT5: ‘At the end of each trial this .rou= 
tine accumulates statistics for inclusion in the run sum- 
mary. If desired the trial summary is displayed. 

Entry INOUT6: At the end of each event statistics 
pertaining to the event are stored and if desired the trial 
status is printed. 

Entry INOUT7 : At the end of each run this routine 
calls subroutine STAT to compute the means and variances of 
each of the data groups compiled during the run, and prints 
the run summary. 

Entry INOUT8: Before each trial this routine sets 
all key variables to their initial values, generates the 
times of the first expected detection events using subrou- 
tines CPAT, DETECT, FIND, and CZDET. The subroutine ELIST 
is called to establish the event calendar for the trial. 

UPDATE: This subprogram computes the attacker and tar- 
get coordinates at the time of the eurrent event using the 
coordinates of the ships at the time of the previous event 


and the X and Y velocity components. With the new 
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positions the range, rate at which range is changing and 

the status of the attacker with respect to the convergence 
zone are calculated. Subroutines DETECT and PDET are 

called to provide the current detection probability estimate. 

FIND: This routine computes the expected counterdetec- 
tion time based on the relative speed vector and counter- 
detection range, R, during trial initialization and after 
each event which results in a course or speed change for 
either unit. (Fig. 6) If present range is greater than R 
and range is increasing counterdetection will not occur. 
When the range is closing from a distance greater than R, 
then CPA range and CPA coordinates are computed. If CPA 
range, Re, is greater than R a detection will not occur. If 
Re is less than R the time of counterdetection is determined 
to be the time at which range between the units will equal 
R. If present range is less than R counterdetection is 
immediate. 

CPAT: This routine computes the time required for the 
attacker to achieve its closest point of approach based on 
the relative speed vector during trial initialization and 
after each event which results in a course or speed change 
for either unit. The time thus computed will be used as 
the time of the next detection probability event. If the 
range is increasing following a course change, but had been 
decreasing, the CPA time is current problem time. If the 
range is decreasing following a course change the coordinates 


of the closest point of approach relative to the attacker 
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Figure 6. Determination of Time Until Counterdetection 


Using Subroutine FIND 
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are determined and the time required for the attacker to 
reach that point is computed based on the relative velocity 
jeter: 

DETECT: This routine evaluates the sonar equation as 
the preliminary step to computing counterdetection time and 
probability. The parameters for the sonar equation are 
compiled using the self and radiated noise curves and the 
ship speeds. The ship noise data is combined with ambient 
noise and sonar gain to form the figure of merit. If a 
detection probability is being computed, the figure of merit 
is returned to the calling program for use with subroutine 
PDET. If counterdetection time is to be established, the 
range of zero signal excess is determined from the propa- 
gation loss curve and is returned to the calling program 
for use in subroutine FIND. 

CZDET: This routine computes the time until the next 
CZ event during each course or speed change by either unit. 
Figure 7 is a logic flow diagram of CZDET. The convergence 
zone is defined in terms of its inner, central and outer 


radii. CZ event times are determined as follows: 


ee Range between the ships is greater than the outer 
CZ radius, and | . 

a. CPA range is less than the central CZ radius: 
the event occurs when the range equals the cen- 
tral CZ radius, 


b. CPA falls between central and outer CZ radii: 


the event occurs when the ships reach CPA, or 
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Figure 7. Subroutine CZDET 
Variable Definitions: 
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Figure 7. (cont. ) 
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c. CPA range is greater than the outer CZ radius: 
no event will be scheduled. 

2, Attacker is within the convergence zone, and 

a. A CZ event has already occured during this 
period within the convergence zone: no event 
will be scheduled, 

b. Range is increasing and greater than the central 
CZ radius: the event time is current problem 
time, : 

c. Range is greater than the central CZ radius and 
CPA range is less than central CZ radius: the 
event will occur when the range equals the cen- 


tral CZ radius, 


Q 


Range is decreasing but CPA range is greater 
than central CZ radius: the event will occur 
when the ships are at CPA, 

e. Range is increasing and less than the central CZ 
radius: the event will occur when the range 
equals the central CZ radius, 

f. Range jig: decreasing and less than the central Ca 
‘radius: the event time is current problem time. 

2 Range is less than the inner CZ radius: the event 
will occur when the attacker reaches the center of 
the convergence zone. 

PDET: This routine computes the probability of detec- 

tion at the beginning of each event, when the attacker 


arrives on search station and when the target changes speed 
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at a new leg time or during evasion, based on the figure 

of merit and range between the ships. The propagation loss, 
L, at the given range is compared with the figure of merit, 
FM, to determine signal excess. Under the assumption that 
propagation loss is a normally distributed random variable 
with known standard deviation, s, then z = (L-~FM)/s repre- 
sents the number of standard deviations by which the propa- 
gation loss exceeds the figure of merit. The probability 
of detection is given by the probability that a standard 
normal random variable is greater than z. The value z is 
compared with the values given in a standard normal cumula- 
tive distribution table to obtain the probability estimate. 

AVOID: This routine computes the course and speed used 
in response to a counterdetection. Figure 8 is a flow 
diagram of subroutine AVOID. Figures 2 and 3 illustrate 
two of the evasion patterns. The calling sequence for a 
special evasion routine is included and may be executed 
through normal inputs. The special evasion must be provided 
by the user as subroutine SPCL. 

INTEL: This routine determines the course, speed and 
position resulting from a surveillance facility intelligence 
detection. Figure 9 shows the logic used in computing the 
intelligence estimate according to the desired option. All 
random error estimates are symetric about the actual value 
with the distribution parameter representing either the 
maximum error in the case of the uniform distribution, or 


two standard deviations of the normal distribution. 
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Figure 8. Subroutine AVOID 
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Figure 9. Subroutine INTEL 
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Computation of course and speed estimates based on position 
error is shown in Figure 1. The first estimate is random 
from a uniform distribution. The first estimate da a new 
leg uses current and previous intelligence positions. Other 
estimates are based on current position estimate and the 
first position estimate for the leg. CZ intelligence is 

not used. 

INTDET: This routine calls INTEL to prepare an 
intelligence estimate, then stores the information received. 
If target course and speed are to be computed based on 
position errors, the base data is generated for the first 
intelligence detection and updated after each target course 
change. If desired intelligence data is printed. 

PROINT: This routine provides the attacker with the 
latest data which is available at a communications period 
but has not been previously transmitted or superceded by a 
CZ intelligence detection. If any such data exists, the 
position estimate is advanced to the current time according 
to the estimated course and speed. Attacker speed is set 
to transit speed, Figure 10 is a flow diagram of subroutine 
PROINT. | 

CZINT: This routine generates the convergence zone 
intelligence data when required. The course, speed, range 
and bearing errors are selected randomly from uniform distri- 
butions using subroutine RAND. Range and bearing errors are 
converted into X and Y coordinates representing the estimated 


target position. 
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Figure 10. Subroutine PROINT 
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ELIST: This routine initializes the event calendar. 
The event calendar consists of a list of event times and a 
list of corresponding event types. The event calendar contains 
all target track pameees all intelligence detection times, 
the times of the next communications event, the next detec- 
tion probability event, the next counterdetection event, the 
next attacker course change event, the next CZ event and 
the mandatory end of trial event. With times assigned to 
all events the lists are arranged in chronological order. 
RPLACE: This routine inserts revised event times in 
the event calendar. At various times throughout the execu- 
tion of a trial Detection, Counterdetection, Attacker 
Course Change, Convergence Zone and eehunt aeons event 
times must be revised. When a revised time is generated, 
the event calendar is searched for the previously scheduled 
time of the event and it is replaced with the revised time. 
The revised time is then sequenced within the event calendar. 
STAT: End of run sample means, M, and sample variances, 
V, are computed for all statistical data compiled during the 
run. Zero valued data may be ignored as in the case of the 
times of ascii detection probability for trials in which 
no positive probability of detection was generated. 


Sample mean is 





where n is the number of data points, x, » being averaged. 
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Sample variance is 


- 2 
) (x, - M) 
k=1 + 


MO n-1l 


CURVE: This is a table look up routine. A table con- 
sists of two lists of corresponding entries. Given a value 
in the range of one of the lists, the corresponding value 
is selected from the other list. When the calling argument 
does not correspond to a table entry, a linear interpola- 
tion between the adjacent entries is performed to produce 
the corresponding value. If the calling argument is outside 
the range of the list, the segment defined by the two data 
points nearest the input argument is assumed to extend 
indefinitely and the output value is determined accordingly. 
Figure 11 illustrates a typical propagation loss curve. 


R R, and R, represent calling arguments with L,, L, and 


1° 2 2 

L, the returned values. In Figure le F, and F, are the 
calling arguments returning Ry and Ro» the expected counter- 
detection ranges corresponding to the given figures of merit. 
Radiated and self noise curves Siways use speed as the 
esis Sesunene: 

ORDER: This routine arranges the event times in chrono- 
logical order. On initializing for each trial and each time 
an event is rescheduled. There are three entry points for 
order. 


Entry ORDER: This section is used to order the en- 


tire list during initialization. For this operation the 


D2 
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Figure 11. Examples of estimating propagation loss relative 
to a specified range using subroutine CURVE. The X's repre- 
sent input data points. R, is less than the minimum input 


range, so L, is determined by extending the first line seg- 


ment. R, is within the convergence zone. 


Propagation 
Loss (db) os 





R Ro —~ Range emacs) 


Figure 12. Examples of estimating detection range based on 
ship noise characteristics. FM, and FM, are figures of 
merit for different ship speeds. The detection range is 


the minimum distance associated with each figure of merit. 
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list is searched for the smallest valued entry. This entry 
and its corresponding event type are exchanged with the 
entries which have index valve 1. The procedure is repeated 
for each index until the ordering is complete. 

Entry ORDER2: This routine places the time of the 
next execution of the event type currently being processed 
in its proper sequence in the event calendar. In this 
case the new event time is located in its proper order and 
the indices of all unexecuted events with earlier times are 
reduced by one. 

Entry ORDER3: This section arranges the time of a 
rescheduled event in its proper sequence in the calendar. 
The event is located in its proper sequence and adjusts 
the intervening times one place up or down as the revised 
time is earlier or later. 

PIM: Attacker course and time to intercept the target 
are computed based on intelligence data each time the 
attacker receives new intelligence, either through a CZ 
event or communications event. The et eleepe time, T, is 


found by simultaneously solving the equations 


So oY = VE 
y 

s,7 + X= vt 
2 2 


where oS and a are the intelligence velocity components, 
Me and Vy are the (unknown) components of attacker velocity, 


‘V, and X and Y are the component differences between the 
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attacker and intelligence positions. If a real, non-nega- 
tive solution exists the intercept coordinates are deter- 
mined and subroutine COURSE is called to find the intercept 
course; 

RAD: This routine converts angles in degrees to radians. 

DEG: This routine converts radians to degrees. 

COURSE: The resultant direction of two component vec-— 
tors is found by this routine using the FORTRAN library 
routine ATAN. 

REVENT: This routine generates random intelligence 
times for each leg of the target track. <A random number, R, 
is selected by the random number generator and compared with 
the probability of an intelligence detection, Pi» associated 
with target track leg i. This comparison is repeated for 
each successive interval of specified duration, t. When R 
is less than Ps the time of detection is given by 
T = T. + tR/P, » where T, is the initial time of the interval. 

RNG: The random number generator selects a random number 
from a uniform distribution between zero and one. This 
routine is related to the shine. Sonb td ion and must be 
rewritten for computers other than the IBM 360. 

NORMAL: This routine is used to generate random numbers 
from a normal distribution using the principal that the sum 
of indepentent random variables tends to a normal distribu- 
tion with mean equal to the sum of the means and variance 


equal to the sum of the variances. Twelve random numbers 


25 





generated by RNG are added to -6 to approximate a normal 
distribution with zero mean and unit variance. 

RAND: A random number, D, is selected uniformly over 
the interval minus one to one using a random number, R, 
generated by RNG according to the formula D = 2(R-0.5). 

SPCL: SPCL is a blank evasion subroutine which may be 
programmed by the user to introduce revised evasion tactics 
or assemble additional evasion statistics. The calling 
sequence is included in the program and controlled by input. 

| FORTRAN library routines: Five subprograms from the 

FORTRAN. library are used by VIM: 

SIN: This routine computes the trigonometric sine of 
an angle given in radians. 

CoS: This routine computes the trigonometric cosine of 
an angle given in radians. 

ATAN: This routine computes the angle in radians whose 
tangent is the calling argument. 

SQRT: This routine computes the square root of a non- 
negative number. | 

TBCOM: This routine controls the input and output func- 


tions of the computer system. 
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IV. SNPUT ff QUEPUD 


A large proportion of VIM is devoted to input of the 
required data, assembly of statistics and information dis- 
play. This chapter is a discussion of the details of the 
input and output routines including the simplified input 
of successive runs and options available for controlling 


the level of output detail. 


A. DATA INPUT 

Input data is treated in six groups: environmental, 
identification, option, situation, target track and intel- 
ligence. This section outlines the input requirements for 
each data group when executing a series of runs. 

Environmental data: Only one set of environmental data 
may be applied to a series of runs. The data must be en- 
tered prior to the first run and remains unchanged through- 
out the series. 

Run identification data: The run identification number 
and description must be included with each Suit Provision 
is made in this group to initiate changes to the option 
and situation data. 

Option aes and situation data: Complete sets of option 
data and situation data must be included preliminary to the 
first run. Either eeaup may be altered between runs, total- 
ly revised or left completely unchanged. Once a change nas 


been entered the new value is retained indefinitely. 
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Target track data: A complete target track must be 
included in the data set for the first run. The target 
track may be totally revised between runs or repeated in 
unaltered form. . 

Intelligence data: Intelligence data may be omitted 
from all runs or included with any of the run data sets. 
Input data may be carried over from run to run as long as 


a random intelligence option is not used. 


Bs -QUTPUT 

Several levels of output are available to provide for 
a careful verification of program operation, to provide 
access to data the analysis of which was not anticipated by 
the design of VIM and to permit efficient computer usage 
when the above factors are not dominant. There are three 
output categories: trial history, trial summary and run 
summary. The trial history and summary may be included with 
the output or deleted. The run summary will always be 
included but may be displayed in an abridged form. 

Trial history: The detailed trial history provides 
a list of each event of the trial with all amplifying data. 
Fach event is characterized by the event time and a brief 
description of the event type followed by the current status 
of detection probability, courses, speeds, range and tacti- 
cal situation along iat any special information which 
might apply to the particular event. The trial history may 
be required for the first trial only, all trials or omitted 


entirely. The trial history adds from one to two seconds 
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of execution time for each trial and adds one output page 
for each five events. 

Trial summary: The trial summary is a list of the 
detection probabilities established during the trial with 
the associated detection timcs and the status of the attack- 
er (transiting or on search station) at each detection 
event. The trial summaries may be displayed with each trial 
history, in place of each trial history, or be omitted 
entirely. When appearing in place of the trial histories 
one page of output is required for each trial. 

Run summary: The run summary is a compilation of the 
individual trial results of some of the significant statis-— 
tics along with ene aiean and variance of each set of statis- 


tics. The individual trial results may be omitted. 


C. DETAILED DESCRIPTION OF OUTPUT 

VIM output includes a display of the input data used 
for the experiment followed by the information generated 
during the run in the form of trial history, trial summary 
and sun summary. 

A ‘Input Data Display (Figure 13) 

All input data is displayed at the beginning of each 
run regardless of the output or input options. The data 
groups shown on the sample input display, Figure 13, are 
explained in the following paragraphs. 

The self and radiated noise curves for the attacker 


and target are each approximated by six linear segments. 
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TRIAL SUM, INPUT INTEL, CZ INTEL 
2 NOPT = 0 NCYCR 
NC = 0O NSUP = 0 NRTT = 
YTI = 500.00 TL = 1.20 
= 0.0 SI = 2.50 
CE = 30.00 SAF = 3.50 
CZCE 15.00 CZBE = 2.50 
SP AX SPAY 
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The first row of each of the noise tables is a list of the 
ship speeds for which noise estimates are to be provided. 
The second row represents the radiated noise in decibels 

at the respective speeds. The third row represents the self 
noise in decibels at the speeds listed in the first row. 

The propagation loss curve is also approximated by 
six linear segments. The top row of the table is the propa- 
gation loss in decibels associated with each of the ranges 
in miles listed in the second row, 

The input value of each option and situation variable 
is listed with the variable name. Each pair of printed lines 
displays the entries from a single input data card. 

The proposed target track is displayed as a sequence 
of rows, each row specifying the data associated with each 
successive leg. The variable names are listed at the top 
of each column. 

Intelligence data input by the user will be dis- 
played as a sequence of rows, each row specifying the data 
associated with a particular intelligence detection. The 
variable names are listed at the cop of each column. If 
intelligence times are determined randomly the times for 
the first trial will appear in rows of seven columns. 
Intelligence times generated randomly for subsequent trials 
will be displayed following.the trial summary. 

2. Detailed Trial History (Figure 14) 
The trial history is sequence of messages that 


describe the state of the trial at each event. The 
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following paragraphs explain the messages found in the 
trial history. 

Event definition message: This message provides 
the problem time of the event and a descriptive phrase iden- 
tifying the event type. The event definition messages are 
self explanatory and refer to the events described in Chap- 
ter III, Section B. Following each event definition is a 
series of status messages which are explained below. 

INTELLIGENCE DATA: The estimate of target position, 
course and speed generated by a surveillance system 
detection. 

CURRENT INTELLIGENCE TIME: The time of the latest 
surveillance intelligence provided to the attacker. If the 
most recent intelligence is the result of a CZ detection the 
time printed is zero. 

EXPECTED CPA TIME: The expected time of CPA based 
on current course and speed is listed if the range between 
the ships is decreasing. 

RANGE OPENING: The range between the units is 
fnereasing 

CURRENT PROBABILITY OF DETECTION: The sonar equa- 
tion estimate of detection probability based on current 
range and speeds. If a speed change by either ship has 
occured during this event the greater weobab lites based on 
the speeds before and after execution of the event is listed. 
The probability shown will not be included in the cumulative 


trial estimate if the attacker is on search station. 
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CUMULATIVE PROBABILITY OF DETECTION: The aggregate 
probability of detection as described in Chapter III, Sec- 
tion A. Note that if the attacker is on search station 
cumulative detection probability will not change until new 
intelligence is received or the trial is terminated at 
which time the maximum probability generated during the 
search period will be used to update the trial estimate. 

ATTACKER RECEIVES LATEST INTELLIGENCE: New intel- 
ligence has become available to the attacker during a Commu- 
nications event or CZ event. 

UPDATED INTELLIGENCE COORDINATES: Intelligence 
data being released to the attacker at a communications 
period is updated to the current time based on the course 
and speed estimates. If the message accompanies a conver- 
gence zone event the values given represent the current 
estimates of a CZ intelligence detection. 

ATTACKER NOT ON STATION: Attacker is awaiting its 
first intelligence information or has not yet reached the 
projected intercept point in response to intelligence data. 

ATTACKER ON STATION: Attacker is on search station 


as determined by the most recently received intelligence. 


ATTACKER COURSE: Self explanatory 
TARGET COURSE: Self explanatory 
TARGET COORDINATES: Self explanatory 


ATTACKER COORDINATES: Self explanatory 
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3. Trial Summary (Figure 15) 

In addition to the cumulative detection probability 
for the trial, each estimated probability with its corre- 
sponding time is listed. The first group of numbers repre- 
sents the times at which a detection probability was 
estimated while the attacker was in transit. JZero's are 
inserted as space keepers. The third group of numbers is 
a list of the times at which a detection probability was 
estimated while the attacker was on search station. Again, 
zero's are inserted as spacers. Sandwiched between the two 
lists are the probability estimates which correspond to the 
times listed above and below. The time of each estimate is 
the non-zero entry in one of the two time lists. 

The random nunber seed, NCYCR, in effect at the end 
of the trial is listed. This number may be used in the 
duplication of the next trial. 

4, Run Summary (Figure 16) 

The run summary lists the data accumulated over a 
series of trials. The following information is provided to 
clarify the messages and descriptive titles appearing in the 
run summary : 

OPTIMAL TARGET INTERCEPT: The time and position 
provided with this message represent the intercept time and 
location that would be generated if the attacker knew the 
target course and speed on each leg and was able to observe 


each course change. 
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Figure 15. Trial Summary 


DATA SET 101 TRIAL 3 PAGE 9 
COMMENCE TRIAL 4 
ATTACKER COURSE = 0.0 SPEED = O.1 
TARGET COURSE = 159.1 SPEED = 14.7 
TARGET COORDINATES —64.37 500.00 
RANGE = 504.13 MILES 
PROBLEM TIME = 74.96 END OF TRIAL 4 


END OF TRIAL SUMMARY OF DETECTION PROBABILITIES - TRIAL 


CUMULATIVE PROBABILITY OF DETECTION 0.147 
TRANSIT PROBABILITY TIMES 
15.64 0.0 
ESTABLISHED PROBABILITIES 
0.07 0.08 
ON STATION PROBABILITY TIMES 
9.0 23.70 


NCYCR = 1452704771 








DATA SET 101 TRIAL 7 PAGE 13 


END OF_GAM 
TOTAL OF 7 «TR 


DATA SET 101 SUMMARY 


INITIAL RANGE TG TARGET 500.00 MILES 


OPTIMAL TARGET INTERCEPT UNDER PERFECT INFORMATION 
TIME 13.89 POSIT -3.090 277.78 


COUNTERDETECTIONS OF ATTACKER BY TARGET EACH TRIAL 
DURING TRANSIT 
0 0 


1 0 0 1 
SAMPLE MEAN = 0.2686 SAMPLE VARIANCE = 0.238 
DURING SEARCH | 
0 0 0 0 0 0 0 
SAMPLE MEAN = 0.0 SAMPLE VARIANCE = 0.0 
TOTAL 
) 0 0 1 0 0 1 
SAMPLE MEAN = 0.286 SAMPLE VARIANCE = 0.238 
pines, te END OA 112.08 74.96 114.60 89.53 75.99 
SAMPLE MEAN = 90.989 SAMPLE VARIANCE = 556.875 
CUMULATIVE PROBABILITIES EACH TRIAL 
is 0.15 0.0 0.15 0.0 0.0 0.15 
SAMPLE MEAN = 0.064 SAMPLE VARIANCE = 0.006 
CUMULATIVE TRANSIT PROBABILITIES EACH TRIAL 
0.0 0.0 0.0 0.07 0.0 0.0 0.03 
SAMPLE MEAN = 0.914 SAMPLE VARIANCE = 0.001 
CUMULATIVE STATION PROBABILITIES EACH TRIAL 
-0 0.15 -0 0. 08 O. 0.0 Osi2 
SAMPLE MEAN = 0.051 SAMPLE VARIANCE = 0.004 


Figure 16. Run Summary 
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DATA SET 101 TRIAL 7 PAGE 14 


r THE FOLLOWING DATA APLIES TO THE TIME JF MAXIMUM OBSERVED PROBABILITY 


NUMBER OF OCCURENCES OF MAXIMUMDETECTION PROBABILITY DURING: 


TRANSIT 0 
SEARCH 
EVASION 0 
MAXIMUM DETECTION PROBABILITIES 
0.0 0.15 0.0 . 0.08 0.0 0.0 0.12 
SAMPLE MEAN = 0.051 SAMPLE VARIANCE = 0.004 
MAXIMUM SEARCH PROBABILITIES 
0.0 0.15 0.0 0.08 0.0 0.0 0.12 
SAMPLE MEAN = G.051 SAMPLE VARIANCE = 0.804 
MAXIMUM TRANSIT PROBABILITIES 
0.0 0.0 0.9 0-07 6.0 0.0 9.93 
SAMPLE MEAN = 0.014 SAMPLE VARIANCE = 0.001 


TIMES OF MAXIMUM DETECTION PROBABILITIES 


0-0 51.29 0.0 23.70 0.9 0-0 23.70 
SAMPLE MEAN = 32.898 SAMPLE VARIANCE = 253.817 
TOTAL POSITIVE ENTRIES 3 
TIMES OF LAST BROADCAST BEFORE MAXPROBABILITIES 
0.0 48. . 2.220 0.0 0.0 2.20 
SAMPLE MEAN = 17.533 SAMPLE VARIANCE = 705.333 


TOTAL POSITIVE ENTRIES 3 


ACQUISITION TIMES OF INTELLIGENCE LEADING TO MAXIMUM PROBABILITIES 


0.0 47-00 0.0 1.00 0.9 Q.0 1.00 
“SAMPLE MEAN = 16.333 SAMPLE VARIANCE = 705.333 
TOTAL POSITIVE ENTRIES 3 
DELAYS FROM ACQUISITION TO BROADCAST 
0.0 1.20 0.0 1.20 0.0 0.0 1.20 
SAMPLE MEAN = 1.200 SAMPLE VARIANCE = -0.990 
TOTAL POSITIVE ENTRIES 3 
DELAYS FROM BROADCAST TO MAXIMUM PROBILITIES 
0.0 3.99 0.0 14.73 9.0 0.0 14.86 
SAMPLE MEAN = 10.895 SAMPLE VARIANCE = 45.648 
TOTAL POSITIVE ENTRIES 3 
CONVERGENCE ZONE DETECTION TIMES 
74.12 0.0 0.0 e 0.0 0.0 0.0 
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SAMPLE MEAN: The arithmetic average and sample 
variance of each category based on the trial results is 
printed under this heading. When the run summary is sup- 
pressed the sample mean and variance only are printed under 
each heading. 

COUNTERDETECTIONS: A list of the number of counter- 
detections each trial, arranged according to the status of 
the attacker, is presented under this title. The first 
list is the number of times the attacker was detected during 
transit, then the number of counterdetections while attacker 
is on search station and, finally, the total. number of 
counterdetections. These lists are omitted if no evasion is 
ordered. 

TIMES TO END TRIAL: A list of the termination time 
of each trial. 

CUMULATIVE PROBABILITIES: A list of the cumulative 
probabilities generated each trial, both overall and accord- 
ing to the attacker status. 

NUMBER OF OCCURENCES OF MAXIMUM DETECTION PROBABI- 
LITY: A count of the number of trials in which the maximum 
aetaerven probability was generated during target evasion, 
attacker transit and attacker search. 

MAXIMUM DETECTION (SEARCH, TRANSIT) PROBABILITIES: 

A list of the maximum detection probabilities achieved each 
trial, both overall and according to attacker status. 

TIMES OF LAST BROADCAST: A list of times of receipt 


of the surveillance intelligence data which preceded the 
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maximum detection probability for each trial. JZero's appear 
for each trial in which the maximum detection probability 
was zero or was preceded by a CZ intelligence detection. 

ACQUISITION TIMES: The times of the surveillance 
intelligence data preceding the maximum detection probabi- 
lity estimates for each trial are listed. If a CZ intelli- 
gence detection preceded the maximum probability, or if the 
maximum probability was zero, a zero entry is printed. 

DELAYS FROM ACQUISITION: The delay from intelli- 
gence acquisition time until broadcast of the final surveil- 
lance intelligence preceding the maximum detection 
probability for each trial is listed. If a CZ intelligence 
detection preceded the maximum probability, or if it was 
zero, a zero entry is printed. 

CONVERGENCE ZONE INTELLIGENCE TIMES: A list of 
times of convergence zone intelligence detections preceding 
the maximum probability of detection. Zero's appear for 
trials in which the maximum detection probability was zero 
or followed a surveillance detection. If CZ intelligence 
is not permitted, this table lists the trial times of CZ 


detection probability estimates which were trial maximums. 


5. Special Messages 


The following messages are printed to indicate 
possible abnormal response to the input data. Each message 
must be interpreted in light of the problem being run. 

UNABLE TO OVERTAKE TARGET DUE TO SPEED DIFFERENTIAL. 


ABORT TRIAL N: The course speed and position estimate upon 
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which the attacker must base its transit path are such that 
an intercept is impossible, and no further surveillance 
intelligence will become available during the trial. This 
message is suppressed “hen the run summary is suppressed. 

INTELLIGENCE TIME OPTION (NOPT) AND INTELLIGENCE 
COMPUTATION OPTION (NCOMP) ARE INCOMPATIBLE. ABORT RUN K: 
An attempt has been made to combine random intelligence 
times with input intelligence estimates. The run with 
identification number K was aborted and the next run exe- 
cuted. The option inputs controlling intelligence should 
be checked. 

SONAR EQUATION PARAMETERS EXCEED INPUT LIMITS. DATA 
EXTRAPOLATION FOR (ATTACKER SPEED, TARGET SPEED, PROPAGA- 
TION LOSS) X: The noise or propagation loss curve was not 
explicitly defined for the value X. A linear extrapolation 
has been performed based on the two data points nearest X. 


This message appears only with the detailed trial history. 
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V. USER'S MANUAL 


The objective of this chapter is to provide a user of 
VIM with all the information necessary to conduct an experi- 
ment with the simulation. Section A describes the inputs 
in detail with cross referencing to pertinent sections of 
the thesis as well as to other input variables for which a 
relationship exists. Tables 2 through 8 provide an input 
guide for the experienced user containing only brief defini- 
tions and no cross referencing. Section C describes the 
assembly of a complete data deck in detail. 

References to inputs in this chapter will. be of the 
form NAME, XX-N, where NAME is the designation of the input 
variable, XX is the input category, and Nis. the serial 
number of the input within the category. The input designa- 


tions and categories are 


EC - Environmental curve, 
RI - Run ddenterreation: 
OD - Option data, 

SD —- Situation data, 

TD - Track data, and 


ID - Intelligence data. 


A. DETAILED DESCRIPTION OF INPUTS 
This section lists each input parameter to VIM by cate- 


gory, and includes the key punching instructions and 
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complete definition. The relationship of each parameter to 
other input variables and its function in the simulation is 
discussed. The column headings found in the input lists 
are number, name, card field, units, limits and description. 
"Number" is the serial number of the input within its category. 
"Name" is the designation of the variable within the program. 
"Card" is the sequence number of the card within the input 
category. "Units" refers to the dimension of the variable: 
miles (mls), hours (hrs), knots kts), decibels (db) or 
degrees (deg). "Limits" defines the range of permissible 
values which may be assigned. If "pos" appears under "limits," 
the variable must be assigned a positive value. If "non-neg" 
is listed, a positive or zero value must be assigned. In 
some cases the limits are determined by the values assigned 
- to other input variables. "Description" includes the var- 
liable definition and amplifying remarks. Columns which do 
not apply to a particular input category are omitted. 
1. Environmental Data Cards 

The environmental data group consists of a title 
card in which the user is free to enter any message for 
reproduction on the input data display, a group of three 
attacker noise cards, a group of three target noise cards 
and a pair of propagation loss cards. Except for the title 
ecard each card is divided into seven fields of ten columns 
each. One entry must be included in each field with a 
decimal point in the appropriate location. Fach entry is. 


indexed according to the position of its field on the card. 


TT 





The entries in the card groups represent point 
approximations of the noise curves. Figures 17 and 18 are 
examples of an attacker radiated noise curve and a propaga~ 
tion loss curve with the input values indicated. Figure 13 
shows the input display of these curves. Interpretation and 
use of the input data points is described in Chapter III, 
Section A under detection, and Chapter TII, Section C under 
DETECT. Interpolation between pairs of data points is 
accomplished by subroutine CURVE. One environmental data 
group must be included with each series of runs. The envi- 
ronmental data is summarized in Table II. 

Title card: The title card is the first card of the 
data set. It may be blank, but must be included. 


Field Description 





2-80 Any title, date or descriptive message. 


Attacker noise curves: 


Name Field Units Description 

VN(1) 1-10 kts Attacker speed card. Each entry is a 

VN(2) 11-20 ' speed for which attacker self and 

VN(3) 21-30 radiated noise must be estimated. Each 

VNC4)-. 31-40 speed listed should be greater than the 

VN(5) 41-50 preceding speed. Entries on the attack- 

VN(6) 51-60 er self noise and radiated noise cards 

VN(7) 61-70 must correspond to the entries on this 
card. 

SNA(1) 1-10 db Attacker self noise card. Entries 

SNA(2) 11-20 represent the self noise generated at 

SNA(3) 21-30 the speeds listed on the attacker speed 

SNA(4) 31-40 card. 


SNA(5) 41-50 
SNA(6) 51-60 
SNA(7) 61-70 
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Figure 17. Example of attacker radiated noise curve approx~ 
4mation. The numbers on the horizontal and 
verticle axes represent the points on which the 


7 approximation is based. See Figures 13 and 20. 
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Figure 18. Example of Propagation loss curve with conver- 
gence zone. The numbers on the horizontal and 
verticle axes represent the points on which the 


approximation is based. See Figures 13 and 20. 
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1--10 
11-20 
21-30 
31-40 
41-50 
51-60 
61-70 


db 


Field Units Description 


Attacker radiated noise card. Entries 
represent the noise radiated by the 
attacker at the speeds listed on the 
attacker speed card. 


Target noise curves: 


Name 








Field Units Description 


1-10 
11-20 
21-30 
31-40 
41-50 
51-60 
61-70 


1-10 
11-20 
21-30 
31-40 
41-50 
51-60 
61-70 


1-10 
11-20 
21-30 
31-40 
41-50 
51-60 
61-70 


kts 


db 


Target speed card. Each entry is a 
speed for which target self and radiated 
noise must be estimated. Each speed 
listed should be greater than the pre- 
ceding speed. Entries on the target 
self noise and radiated noise cards must 
correspond to the entries on this card. 


Tareet self noise card. Entries repre- 
sent the self noise generated at the 
speeds listed on the target speed card. 


Target radiated noise card. Entries 
represent the noise radiated by the tar- 
get at the speeds listed on the target 
speed card. 


Propagation loss curve: 


Field Units Description 


1-10 
11-20 
21-30 
31-40 
41-50 
51-60 
61-70 


ab 


Propagation loss card. Entries repre- 
sent the reduction in noise level at the 
distances from the noise source listed 
on the range card. No restrictions are 
placed on the relative size of the 
entries, but decreasing entries associ- 
ated with a convergence zone will result 
in jump increases in counterdetection 
range at the point of decrease. See 
Figure le. 
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Name Field Units Description 





RDET(1) 1-10 mls Range curve. Entries represent the 


RDET(2) 11-20 distance from the noise source associ- 
RDET(3) 21-30 ated with each entry on the propagation 
RDET(4) 31-40 less card. Entries should be 

RDET(5) 41-50 increasing. 


RDET(6) 51-60 
RDET(7) 61-70 


2. Run Identification Cards 
The run identification cards must be included with 
each run of a series of runs. The first card is the run 
number card with one to four integer entries which must be 
right adjusted in their fields, and contain no decimals. / 
The second card is the run description card. The run iden- 
tification cards are summarized in Table III. 


Run number card: 


No. Name Field Description 


1 NUMBER 1-5 An arbitrary run identification number 
assigned by the user. NUMBER is printed 
with the input data display, the run sum- 
mary and at the top of each page of output. 


2 NCYCR 6-15 Random number seed. This entry is blank 
“on the first run of a series of runs, and 
whenever a new option data card is includ- 
ed with the run data. In that event, NOC, 
RI~3,is zero or blank and NCYCR is entered 
as OD-6. 


3 Nc 16-20 Number of option changes this run. This 
entry is zero or blank on the first run 
and whenever a new option data card is 
included with the run data. If NOC has 
the value one through 13, option change 
ecards specifying NOC changes to the option 
data are required. If no change to the 
option data is desired set NOC to 99 and 
omit option data and option change cards. 
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No. Name Field Description 


4 NC 21-25 Number of situation changes this run. 
This entry must be zero or blank on the 
first run and whenever a new set of situa- 
tion data is to be entered. If NC has 
the value one through 27, situation change 
cards specifying NC changes to the situa- 
tion data are required. If no change tO 
the situation data is desired set NC to 
99 and omit situation data and situation 
change cards. When the option data card 
is to be included with the data set NC 
is entered as OD-10. 


Run description card: The run description card 
contains alphabetic and numeric information or may be blank 


but must be included with the data set for each run. 


Field Description 
2-80 A date, title or descriptive message 


regarding the run. 


3. Option Data Cards 


The option data card must accompany the first run 
data set and may be omitted, replaced by option change cards 
or included with subsequent data sets. All inputs are inte- 
ger and must be right adjusted within the specified field. 
An option input card must be included whenever NOC, RI-3, 
is zero or blank. The option data inputs are summarized in 


Table IV. 


No. Name Field Limits Description 


1 NTRIAL 1-5 1-100. Number of trials for this run. 
Chapter VI, Section A includes a 
discussion of determining the mag- 
nitude of NTRIAL. 
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No. Name Field Limits Description 


2 M ° 6-10 
3 N 11-15 
i TORE 16-20 


5 NOPT 21-25 


0-50 


0-1.00 


1-5 


HOST 


The number of target track cards 

included with the data set for this 
run. If M is zero the target track 
from the previous run will be used. 


The number of intelligence detec- 
tions to be input for this run. If 
N is zero intelligence detection 
times will be taken from the previ- 
ous run or computed randomly per 
NOPT, OD-5. 


Level of output detail. See Chap- 
ter V, Sections B and C for fur- 
ther discussion. The run summary 
is included with all output 
options. 


TOPT=1: Detailed history of each 
trial including the end of trial 
summary. If random intelligence 
times are computed for each trial 
(NOPT, OD-5) they are printed at 
the start of each trial.. 


IOPT=2: Detailed history of first 
trial. End of trial summary for 
each trial. If random intelligence 
times are computed for each trial 
(NOPT, OD-5) they are printed pre- 
ceding the trial summary display. 


TOPT=3: End of trial summary for 
each trial. If random intelligence 
times are computed before each trial 
(NOPT, OD-5) they are printed pre- 
ceding the trial summary display. 


IOPT=4: Detailed history of the 
first trial only. 


IoPpT=5: End of run summary only. 


Intelligence time option switch. 
Table I summarizes the values of 
NOPT, NCOMP (OD-8) and INRAND 
(OD-9) for each of the permissible 
intelligence combinations. 
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No. Name 


Field Limits 


NOPT (cont.) 


6 NCYCR 26-35 1-10" 


Description. 


NopT=-1: Intelligence detection 
times will be selected randomly 
before each trial according to the 
parameters TINT, SD-6, and PINT, 


‘mp4, Under this option N, OD-3, 


must be zero and NCOMP, OD-8, must 
be one or two. Omit the intelli- 
gence data cards. 


NoPT=0: Intelligence detection times 
will be taken from the intelligence 
data cards. Set N, OD-3, to the 
number of intelligence times desired 
and include N intelligence data 
cards. 


NOPT=1: Intelligence detection times 
will be selected randomly for the 
first trial according to the para~ 
meters TINT, SD-6, and PINT, TD-4. 
Subsequent trials will use the same 
intelligence times. Under this 
option N, OD-3, must be zero and 
NCOMP, OD-8, must be one or two. 

Omit the intelligence data cards. 


Random number seed. NCYCR is the 
first random number used by subrou- 
tine RNG in generating a sequence of 
random numbers. RNG produces a new 
value of NCYCR each time it is called. 
When the option data card is omitted 
NCYCR is entered as Ri-e. 


Evasion option. The evasion pattern 
to be followed by the target after 
counterdetection is designated by 
this input. See Chapter III, sec- 
tion A for further discussion. 


NEVD=1: Target reverses course and 
reduces speed by one half. 


NEVD=2: Target alters course clock- 
wise from the true bearing of the * 
attacker, and changes speed. The 
magnitude of the course change is 
defined by EVANG, SD-15. The speed 
change is defined by EVSPD, SD-16. 
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No. Name 


8 NCOMP 


Field Limits Description 


NEVD (cont.) 


46-50 


ae 


NEVD=3: Target turns away from the 
attacker track and steers a course 
at an angle EVANG, SD-15, relative 
to the true bearing of the attacker, 
and changes speed by the amount 
EVSPD, SD-16. Evasion option 3 is 
pictured in Figure 2 


NEVD=4: Target turns away from the 
target track and steers a course at 
an angle EVANG, SD-15, relative to 
the reciprocal of attacker course, 
and changes speed by the amount 
EVSPD, SD-16. Evasion option 4 is 
pictured in Figure 3. 


NEVD=5: No evasion. Counterdetec- 
tions are suppressed. 


NEVD=6: Special evasion. Subrou- 
tine SPCL is called. SPCL is a 

dummy routine which may be programmed 
to meet any special evasion require- 
ments. In its current state, this 
option will result in the recording 
of counterdetections at the minimum 
interval of ETIME, SD-7. 


Intelligence data option switch. 
This option designates the method of 
generation of surveillance intelli- 
gence estimates. Table I summarizes 
the values of NCOMP, NOPT (OD-5) and 
INRAND (OD~9) for each of the per- 
missible intelligence combinations. 


NCOMP=1: Target estimated position, 
course and speed are selected random- 
ly from the distribution designated 
by INRAND, OD-9. The parameters of 
the position estimate are SPAX, TD-5, 
and SPAY, TD-6. The parameter for 
the course estimate is CE, SD-18. 

The parameter for the speed estimate 
is SE, SD-17, with the additive 
factor SAF, SD-19. 


NCOMP=2: Target estimated position 
is selected randomly from the distri- 
bution designated by INRAND, OD-9, 
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No. Name Field Limits Description 


NCOMP (cont.) NCOMP=2 (cont.): 
based on the parameters SPAX, TD-5, 
and SPAY, TD-6. The course and 
speed estimate is based on two posi- 
. tion estimates and the intervening 
° time. Parameters for the initial 
course and speed estimates are CE, 
SD-18, and SE SD-17 along with the 
additive factor SAF, SD-19. Speed 
estimates are constrained to fall 
between S-SE and S+tSE+SAF, where 5 
is the target leg speed. The compu- 
tation process is described in Chap~ 
ter III, Section C under subroutine 
‘INTEL, and illustrated in Figure l. 


NCOMP=3: Target position is derived 
from input position error informa- 
tion taken from the intelligence 
Gata cards. Course and speed are 
selected randomly as described for 
NCOMP=1. Estimated target position 
is on bearing BINP, ID-3, at range 
RINP, ID-2, from the actual target 
position at the time of the intelli- 
gence detection. 


NCOMP=4: Target position, course 

and speed estimates are taken from 
the input information on the intelli- 
gence data cards. Determination of 
the position estimate is as described 
for NCOMP=3. Course and speed esti- 
mates are taken directly from CINP, 
ID-4, and SINP, ID-5. 


NCOMP=5: Target position estimate 

is derived from the position error 
information input on the intelligence 
data cards. The course and speed 
estimates are based on two position 
estimates and the intervening time 

as described for NCOMP=2. Computa- 
tion of the position estimate is as 
described for NCOMP=3. 


Note: If NCOMP is assigned the value three, four 
or five NOPT, OD-5, must be zero. 
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No. Name Field Limits Description 


9 INRAND 46-50 1-4 Intelligence error distribution 
- switch. Table I summarizes the 
values of NCOMP (OD-8), NOPT (OD-5) 
and INRAND for each of the permissi- 
ble intelligence combinations. If 
NCOMP is four or five INRAND is 
ignored.  - 


INRAND=1: Intelligence position 
course and speed errors are selected 
from a uniform distribution. Maxi- 
mum position errors are specified by 
SPAX, TD+5, and SPAY, TD-6, Maximum 
speed error is SE, SD-17, Maximum 
course error is CE, SD-18. 


INRAND=2: Intelligence position 
error is selected from a normal dis- 
tribution. Course and speed errors 
are from a uniform distribution as 
for INRAND=1. The standard devia- 
tions for the * and position 
errors are given by 1/2 SPAX, TD-5, 
and 1/2 SPAY, TD-6. 


INRAND=3: Intelligence position 
error is selected from a uniform 
distribution as described for 
INRAND=1. Course and speed errors 
are selected from a normal distribu- 
tion with standard deviations of 

1/2 CE, SD-18, and 1/2 SE, SD-17, 
respectively. 


INRAND=4: Intelligence position 
error is selected from a normal 
distribution as for INRAND=e2. 

Course and speed errors are from a 
normal distribution as for INRAND=3. 


10 NC 57-55 0-27, The number of changes to situation 
99 inputs from the previous trial. If 

NC is zero a complete set of situa- 
tion inputs is required. If NC is 
99 the situation data from the pre- 
vious run will be used, and the 
situation data and situation change 
cards must be omitted from the data 
deck. If NC is assigned a value 
from one through 27 situation change 
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No. 


Lid 


12 


BRS. 


Field Limits Description 


Name 

NC (cont.) 
NSUP 56-60 
NRTT 60-65 
IcZ 66-70 


Oi: 


Ov 


0,1 


ecards specifying NC changes must be 
included. If the Option card is 
omitted NC is entered as RI-4. 


Run summary supression switch. 
When activated NSUP reduces the 
volume of output displayed in the 
run summary. 


NSUP=0: The full end of run summary 
will be printed. 


NSUP=1: The lists of individual 
trial results will be deleted from 
the run summary. Sample means and 
variances will be retained. 


Random target track option. When 
activated NRTT causes the target 

to deviate randomly from its assigned 
track (target track data). 


NRTT=0: Target track data is used 
as given. 


NRTT=1: Target initial position and 
the course and speed for each leg 
are selected from uniform distribu- 
tions. Deviation of initial target 
X and Y coordinates from the input 
initial position is limited by 
SPAX(1), TD-4, and SPAY(1),, TD=5. 
The input initial position is given 
by XTL, -SD=35, and, TTL, SD-4. Maxi- 
mum course and speed deviations are 
defined by CE, SD-18, and SE, SD-17. 
The base course and speed for each 
leg is CST, TD-1, and SPT, TD-e. 


Convergence zone intelligence option. 
This switch is in effect any time 
CZR, SD-22, is greater than zero. 


ICZ=0: No convergence zone intelli- 
gence will be generated. Detection 
probability estimates generated with- 
in the convergence zone will be 

used in establishing the trial proba- 
bility estimate. 
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No. Name Field Limits Description 


10Z. -Ceonts.) TCZ=1: Probability estimates gene- 
rated within the convergence zone 
are converted to intelligence esti- 
mates which are provided immediately 
to the attacker. If an intelligence 
estimate is to he generated speed, 
“course range and bearing errors are 
selected from uniform distributions 
limited by CZSE, SD-24, CZCE, SD-25, 
CZRE, SD-26, and CZBE, SD-27, 
respectively. 


Option change cards: When NOC, RI-3, is assigned a 
value from one eneough 13 the option data card is replaced 
by one or more option change cards naming NOC changes to 
the option data from the previous run. Up to eight changes 
may be entered on an option change card. Each change is 
entered as a pair of numbers: the first is the option input 
number; the second is the revised value. All entries are 
integer and must be right adjusted within their fields. If 


NOC is zero, blank or 99 the option change card is omitted. 


Field Description 

x1-x5 Option input serial number. 

x6—yO Revised option input value. 

(x represents any digit zero through seven. y is 


x plus one.) 


89 





Table I. Intelligence Option Summary. 

This table summarizes the input values for the 
intelligence option switches NOPT, NCOMP and INRAND to 
achieve each of the possible methods for generating intelli- 
gence estimates. The effect of NCOMP values 1 and 2 is 
independent of the NOPT setting. When NCOMP is assigned 3, 


} or 5, NOPT must be 0. 


NOPT NCOMP INRAND Method of estimation 


=a Random intelligence times each trial. 
0 Input intelligence times. 
1 Random intelligence times first trial. 


1 1 Random uniform position estimate. 
Random uniform course and speed est. 


2 Random normal position estimate. 
Random uniform course and speed est. 


3 Random uniform position estimate. 
Random normal course and speed est. 


yf Random normal position estimate. 
Random normal course and speed est. 
2 m3 Random uniform position estimate. 
Course and speed based on positions. 
2,4 Random normal position estimate. 
Course and speed based on positions. 
0 Input intelligence times. 


3 ali Ae Input position estimate. 
Random uniform course and speed est. 


3,4 Input position estimate. 
Random normal course and speed est. 
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Table I. (cont.) 
NOPT NCOMP INRAND Method of estimation 
0 Input intelligence times. 


y 1-4 Input position estimate. 
Input course and speed est. 


Sa) 
ft 
I 
= 


Input position estimate. 
Course and speed based on 
positions. 


91 





4. Situation Data Cards 
The situation data cards must accompany the first 
run data set and may be omitted, replaced by situation 
change cards or included with subsequent data sets. All 
inputs are decimal and may be placed anywhere within their 
fields. A decimal point must be included with each entry. 
The situation input cards must be included whenever NC, 


RI-4 or OD-10, is zero or blank. The situation data is sum- 


marized in Table V. 


No. Name Card Field Units Limits Description 
1 TEND 1 1-10 hrs pos Nominal time to end the 

trial. No random intelli-~ 
gence times will be generat- 
ed after TEND. The trial 
will terminate at TEND + 100 
if the normal criteria have 
not been satisfied. See 
Chapter III, Section B for 
details. 


2 cINT il 11-20 hrs O-TEND Communication interval. 
Attacker monitors the intel- 
ligence broadeast every CINT 
hours when not on search 
station. While on station 
the attacker monitors the 
broadcast every two hours. 
If CINT is zero the attacker 
will receive intelligence 
whenever it becomes avail- 
able (TL, SD-5). 


3 XTI 1 £21-30 mis Initial X and Y coordinates 

yo YPT ak 31-40 mils of the target relative to 
the attacker. xXTI is the 
east-west displacement. 
YTI is the north~south dis- 
placement. Positive entries 
are north or east. When the 
random target track option 
is selected (NRTT, OD~12), 
(XTI,YTI) represents the 
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No. Name Card Field Units Limits Description 


NOs. NCIS 2 On Sa 


YTI (cont.) center of a rectangle with- 
in which the target is ran- 
Gomly located. 


5 TL 1 41-50 hrs pos Fixed intelligence delay. 
Surveillance intelligence 
is held for TL hours before 
being released for communi 
cation to the attacker. 


6 TINT 1 51-60 hrs pos Intelligence detection in- 
terval. Intelligence detec- 

tions occur at most one 
time during each successive 
time period of duration 
TINT whenever random intel- 
ligence times are specificd 
(NOPT, OD-5). The probabi- 
lity of a surveillance de- 
tection within the interval 
4s PINT, TD-5. The exact 
time of detection is placed 
randomly within the interval. 


7 |EVIME 1 61-70 hrs pos Evasion time. Upon detect- 
ing the attacker the target 


will execute the assigned 
evasion maneuver (NEVD, 
OD-7) for an interval of 
length ETIME. Counterdetec-— 
tions are suppressed during 
the evasion interval. 


8 ss 2, 1-10 kts pos Search speed. The attacker 
uses this speed while on 
station for improved sonar 
performance. 


9 SA 2 11-20. kts pos Transit speed. The attacker 
uses this speed to proceed 
to search station in res- 
ponse to intelligence data. 


10 VI 2 21-30 kts pos Initial speed of attacker 
used prior to receipt of 
the first intelligence data. 


12 -CALE 2 31-40 deg 0-360 Initial course of attacker 
used prior to receipt of 
the first intelligence data. 
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No. Name Card Field Units Limits Description 
a2. Str 2. 41-50 hrs pos Search leg time. Upon 

arrival on station the at- 
tacker searches on a course 
perpendicular to the esti- 
mated target course for 

‘i 1/2 SI, then reverses 
course every SI hours there- 
after until new intelli- 
gence is received. 


13 AMB 2 51-60 db Ambient noise level. AMB 
is a term of the sonar equa- 
tion representing a lower 
bound on self noise regard- 
less of the values derived 
from the self noise curves. 


14 SIGMA 2 61-70 db pos Standard deviation of pro- 
pagation loss. SIGMA is 
used with the sonar equation 
to estimate detection 
probabilities. 


15 EVANG 3 1-10 deg 0-360 Target evasion angle. In 
executing evasion patterns 
2, 3 or 4 target evasion 
course is determined by the 
geometry at counterdetection 
and EVANG. (NEVD, OD-7) 


16 EVSPD 3 11-20 kts note Evasion speed increment. 
In executing evasion pat- 
terns 2, 3 or 4 target eva- 
sion speed is leg speed 
plus EVSPD. (NEVD, OD-7) 


Note: EVSPD may be positive or negative, but if ne- 

gative it must not be permitted to drive the 
target speed to zero or below during evasion. 
Target leg speed, SPT, TD-2, random track op- 
tion NRTT, OD-12, and speed error, SE, SD-17, 
must be considered. 


17 SE 3 21-30 kts note Speed error. Intelligence 
speed estimation by the sur- 
veillance facility uses SE 
as the maximum error of a 
uniform distribution or two 
standard deviations of a 
normal distribution. In 
the latter case SE defines 
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Sasa CREED AnIOCRNC Ca 


SE (cont. ) a region in which 98.85% of 

. the errors are expected. 

When the intelligence speed 
estimate is computed from 
position estimates SE is 
the maximum positive speed 
error, and is combined with 
SAF, SD-19, to form a limit 
on negative error. SE is 
the maximum deviation from 
assigned target track 

speed SPT, TD-2, under the 
random target track option, 
NRTT, OD-12. 

Note: Under the fixed target track option SE must 
be positive and less than the minimum target 
leg speed. Under the random target track op- 
tion SE must be less than one half the mini- 
mum target leg speed. This will assure that 
target speed and estimated speed will not be 
driven to zero or below. 


18 CE 3 31-40 deg 0-360 Course error. Intelligence 
course estimation by the 
surveillance facility uses 
CE as the maximum error of 
a uniform distribution or 
two standard deviations of 
a normal distribution. In 
the latter case CE defines 
a region in which 98.85% of 
the errors are expected. 

CE is the maximum deviation 
from the assigned target 
track course, CST, TD-3, 
under the random target 
track option, NRTT, OD-le. 


19 SAF 3 41-50 kts pos Safety factor. When the 
surveillance intelligence 
speed estimate is generated 
by VIM, and for all CZ in- 
telligence estimates, SAF 
is added to the estimated 
speed in an attempt to force 
the attacker to a search 
station ahead of the actual 
target. When intelligence 
speed is based on position 
estimates. SAF is combined 
with SE, SD-17, to 
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No. Name Card Field Units Limits Description 
SAF (cont.) establish a limit on nega- 
tive speed error. 


20 GAINA 3 51-60 db Attacker sonar gain. GAINA 

is used in the sonar equa- 
= tion for estimating detec- 

tion probability, and is 
composed of recognition 
differential, Nrd, and di- 
rectivity index, Ndi, asso- 
ciated with the attacker's 
sonar suite) (GAINA=Nrd-Ndi) 


21 GAINT 3 61-70 db Target sonar gain. GAINT 
is used in the sonar equa- 
tion for determining coun- 
terdetection range, and is 
composed of recognition 
differential, Nrd, and di- 
rectivity index, Ndi, asso- 
ciated with the target's 
sonar suite. (GAINT=Nrd-Ndi) 


22 OZR y 1-10 mls non-neg Convergence zone central 
radius. CZ events occur 
when the range between the 
ships is as near the cen- 
tral CZ radius as possible 
based on the current rela-~ 
tive track. If no conver- 
gence zone is represented 
on EC-5,CZR should be set 
to zero. The nature of the 
CZ event is designated by 
ICZ, OD-13. 


23 «CZW 4 11-20 mls 0O-CZR Convergence zone half width. 
The outer CZ radius is 

given by CZW + CZR (SD-22). 
The inner CZ radius is CZR- 
CZW. During the time the 
range between the units 
falls within the minimum 
and maximum CZ radii exactly 
one CZ event will occur. 


24 cCZSE 4 21-30 kts note Convergence zone speed 
error. When a convergence 
zone intelligence estimate 
is generated, the speed 
error is selected randomly 
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No. Name 


Card Field Units Limits Description 


from a uniform distribution 
with CZSE representing the 
maximum error. (ICZ, OD-13) 


To prevent a negative speed estimate, CZSE 


should be less than the minimum target leg 


If the random target track option 


(NRTT, OD~12) is activated the minimum target 
leg speed should be greater than CZSEtSE, 


CZSE (Ccont.) 
Note: 
speed. 
SD-17. 
25 CZCE 4 31-40 deg 0-~360 
26 CZBE 4 41-50 deg 0-360 
27 CZRE 4 


51-60 mis 


Situation change cards: 


Convergence zone course 
error. When a convergence 
zone intelligence estimate 
is generated, the course 
error is selected randomly 
from a uniform distribution 
with CZCE representing the 
maximum error. (ICZ, OD-13) 


Convergence zone bearing 
error. When a convergerice 
zone intelligence estimate 
is generated, the position 
error is computed from a 
bearing and range error 
(CZRE, SD-27). The bearing 
error is selected randomly 
from a uniform distribution 
with CZBE representing the 
maximum error. (ICZ, OD-13) 


Convergence Zone range 
error. When a convergence 
zone intelligence estimate 
is generated, the position 
error is computed from a 
range and bearing error 
(CZBE, SD-26). The range 
error is selected randomly 
from a uniform distribution 
with CZRE representing the 
maximum error. (ICZ, OD-13) 


When NC, RI-4 or OD-10, is 


assigned a value from one through 27 the situation data cards 


are replaced by one or more situation change cards naming NC 


changes to the situation data from the previous run. 
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Up to 





five changes may be entered on each change card. Each 
change is represented by a pair of numbers: the first is 
the integer valued situation input number; the second is 
the revised (decimal) value which may appear anywhere in 
its assigned field and must include a decimal point. If 
NC is zero, blank or 99 the situation change cards must be 


omitted. 


Fields Description 
be oD 
16-20 
31-35 Situation input serial numbers. 
46-50 
61-65 
6-15 
21-30 
36-45 Revised situation input values. 
51-60 
66-75 


5. Target Track Cards 

Course, speed and surveillance system detection 
parameters for each leg are the components of the target 
track data. One card is required for each leg of the tar- 
get's transit path, with the number of legs specified by M, 
OD-2. If Mis zero the target track cards must be omitted 
and the target track from the previous run will be used. 
Target track cards must be included with the first run data 
set. All entries are decimal and may be located anywhere 
within their fields. Decimal points must be included. The 


target track data is summarized in Table VI. 
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No. Name Field Unite Limits Description 


1 CST(I) 1-10 deg 0-360 Base course for leg I. Prior 
. to the leg termination time, 
TT, TD-3, the target will de- 
viate from this course only 
in response to a counterdetec- 
tion (NEVD, OD-7) or in res- 
ponse to the random track 
option, NRTT, OD-le. 


2 SPT(I) 11-20 kts pos Base speed for leg I. Prior 
to the leg termination time, 
TT, TD-3, the target will 
deviate from this speed only 
in response to a counterdetec-— 
tion (NEVD, OD-~7) or in res- 
ponse to the random track 
option, NRTT, OD-12. 


2 PCL) 21230 hrs pos Terminal time for leg I. When 
problem time equals terminal 
time for the leg, the target 
track index is incremented and 
the parameters for the next 
track leg are activated. The 
terminal time for the final 
leg must exceed TEND, SD-1, 
by more than 100 hours. 


Hh =PINT(I) 31-40 0-1 Intelligence detection proba- 
bility. Under the random 
intelligence time option, 
NOPT, OD-5, the surveillance 
facility will detect the tar- 
get on leg I with probability 
PINT during each successive 
interval of length TINT, SD-6. 


SPAX(I) 41-50 mls Surveillance position error 

SPAY(I) 51-60 mls parameters. When a random 
position estimate is generated 
by VIM, SPAX and SPAY are the 
maximum X and Y deviation of 
a uniform error, or two stan- 
dard deviations of a normally 
distributed error. In the 
latter case, SPAX and SPAY de- 
fine regions of the X and Y 
axes in which 98.85% of the 
error is expected. (NCOMP, 
OD-8; INRAND, OD-9) When the 
random track option is 


OVI 
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No. Name Field Units Limits Description 
SPAX, SPAY (cont.) activated (NRTIT, OD-12) 
SPAX(1) and SPAY(1) represent 


the maximum deviations from 
XTI, SD~3, and YTI, SD-4., 


6. Intelligence Data Cards 
The intelligence data cards refer to surveillance 

system detections and consist of detection time and related 
course, speed and position estimates. One card is required 
for each intelligence detection. The number of detections 
is specified by N, OD-3, which should be zero whenever NoPT, 
OD~5, is not zero. If both N and NOPT are Zero, the intel- 
ligence data cards must be omitted and the intelligence data 
input for the previous run will be used. If N is wero and 
NOPT is not zero the intelligence data cards must be omitted 
and random intelligence data will be generated within VIM. 
All entries are decimal and may be located anywhere within 


their fields. Decimal points must be included. The intel- 


ligence data is summarized in Table VII. 


No. Name Field Units Limits Description 


1 TI(J) 1-10 hrs pos The time of intelligence detec- 

. tion J. When problem time 
equals intelligence detection 
time, the course speed and 
position will be estimated in 
accordance with the intelli- 
gence option switches NCOMP, 
OD-8, and INRAND, OD-9. 


2 RINP(J) 11-20 mis_ pos Range error for intelligence 
detection J. If NCOMP, OD-8, 
is one or two this entry may 
be blank. Otherwise, intelli- 

- gence position is based on 
range error and bearing from 
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No. Name Field Units Limits Description 


x RENE Veontg: ~ actual target position (BINP, 
: ID-2). 


3 BINP(J) 21-30 deg 0-360 Bearing from target of intel- 
ligence position estimate J. 
If NCOMP, OD~8, is one or 
two this entry may be blank. 
Otherwise, intelligence posi- 
tion is based on the bearing 
from the target and the input 
range error BINP, ID-e. 


4 CINP(J) 31-40 deg 0-360 Course estimate for intelli- 
gence detection J. This entry 
is required only when NCOMP, 
OD-8, is assigned the value 
four. 


5 SINP(J) 41-50 kts pos Speed estimate for intelli- 
- gence detection J. This entry 
is required only when NCOMP, 
CD-8 is assigned the value 
four. 


7. Terminal Card (99999 Card) 
Following the last data card from the data set for 
the final run, the terminal card signais the end of the data 


set. This card is described in Table VIII. 


Field Description 
1-5 99999 
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Table II. Environmental Data Summary 
Title card: 
Field Description 


2-80 Users message. 
Attacker noise curves: Decimal entries, 


Name Field Units Description 
VN(I) 10 col kts Attacker speed card. I=1,...,7. 
SNA(I) 10 col db Attacker self noise. Be Eee ener ars 


RNA(I) 10 col db Attacker radiated noise. Tayo Ts 


Target noise curves: Decimal entries. 
Name FPiéic Units Description 
STN(I) 10 col kts Target speed card. TA oaks ty 
SNT(I) 10 col ab Target self noise. De We ein Pe 


RNT(I) 10 col ab Target radiated noise. se ae een! ee 


Propagation loss curve: Decimal entries. 
Name Field Units Deseription 
PLOS(T) 10 col db Propagation loss. ig nearer ere Ge 


RDET(I) 10 col mls Range of loss. ee Des pateras 6% 
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Table IIT. Run Identification Summary 
Run number card: Integer entries. 
No. Name Field Description 
1 NUMBER 1- 5 Run identification number. 
2 NCYCR 6-10 Random number seed (OD-6). 
3 NOC 11-15 Number of option changes. 


4 NC 16-20 Number of situation changes. (OD-10). 


Run description card: User's message in columns two 


through 80. 
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Table IV. 
No. Name 

1 NTRIAL 

a 6M 

3 4#ON 

4 JOPT 

5 NOPT 

6 NCYCR 


Option Data Summary: Integer Entries. 


Field 
1-100 
6-10 


11-15 


16-20 


21-25 


26~35 


Limits 
1-100 
O- 50 


0-100 


im 5 


=1,0,1, 


ig?” 


Description 
Number of trials. 
Number of target tracks cards. 


Number of intelligence detection 
ecards. 


Output option. 
IOPT=1: Details each trial; 
trial sum. 
e: Details first trial; 
trial sum. 
3: Trial summaries. 
4: Details first trial. 
5: No trial details. 


Intelligence time option. 
NOPT=-1: Random times each trial. 
O: Input times. 
1: Random times first trial 
for all trials. 


Random number seed. (RI-~2) 


Evasion option. 
NEVD=1: Course reversal. 
2: Course change relative to 
attacker bearing. 
3: Course change away from 
attacker bearing. 
N: Course change away from 
attacker course. 
5: No evasion. 
6: User's choise (SPCL). 


Intelligence data option. 
NCOMP=1: Random position, course 
and speed. 
2: Random position, computed 
course and speed. 
3: Input position, random 
course and speed. 
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11 


le 
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Name Field Limits Description 





NCOMP (cont. ) 4: Input position, course 
and speed. 
5: Input position, computed 
course and speed. 


TNRAND 46--50 Tuell Intelligence error distribution 
switch. 
INRAND=1: Uniform position, course 
and speed. 
2: Normal position, uniform 
course and speed. 
3: Uniform position, normal 
course and speed. 
4: Normal position, normal 
' course and speed. 


NC 51-55 0-27, Number of situation data changes. 
99 (RI-4) 
NSUP 56-60 0,1 Run summary supression switch. 


NSUP 1: Individual trial statistics 
omitted from run summary. 


NRIT 61-65 0,1 Random target track option switch. 
NRTT 1: This option is activated. 
TEZ 66-70 0,1 Convergence zone intelligence 
option. 


ICZ 0: Detection probabilities 
included overall estimate. 


1: Detection probabilities 
converted to Ca 
intelligence. 

Option change card: Integer entries 
Field Description 


x1-x5 Option input number. 
x6-y0 Option input value. : 


(x represents any digit zero through seven. y is x 


plus one.) 
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Table V. Situation Data Summary: Decimal Entries. 


No. Name Card Pield Units Limits Description 


1 TEND 1 1~10 hrs pos Nominal time to end 
trial. 


2 CINT 1 11-20 hrs O-TEND Communications interval. 
CINT O: Intelligence 
transmitted when 
ready. 


ATE 
Vil 


21-30 mls Initial target coordi- 
31-40 mls nates. 


Ww 
HH 


iy 1 47-50 hrs pos Intelligence delay time. 


2 
6 TINT 1 51-60 hrs pos Intelligence detection 
interval. 


7 |ETIMVE 1 61-70 hrs pos Evasion time. 


a -8s 2 1-10 kts pos Attacker search speed. 
9 SA 2 11-20 kts pos Attacker transit speed. 
10 VI 2 21+30 kts pos Initial attacker speed. 


dk SORE 2 31-40 deg 0-360 Initial attacker course. 


12) - 187 2 41-50 hrs pos Attacker search leg 
time. 

13 AMB 2 51-60 db Ambient noise level. 

14 SIGMA 2 61-70 db pos Standard deviation of 


propagation loss. 


3 1-10 deg 0-360 Target evasion angle. 
16 EVSPD 3 11-20 ‘kts note” Target evasion speed. 
17 SE 3 21-30 kts note Intelligence speed error. 
18 CE 3 31-40 deg 0-360 Intelligence course 
error. 
1 


These entries must be chosen so that target track 
speeds or intelligence speed estimates cannot take on nega- 
tive values. 
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No. Name Card Field Units Limits Description 


: 19 SAF 3 41-50 kts pos Intelligence speed 
safety factor. 

20 BAINA 3 51-60 ° db Attacker sonar gain. 
21 GAINT 3 61-70 db . Target sonar gain. 
22 °CZR y 1-10 mls non-neg CZ central radius. 
23 OC AW 4 11-20 mls O-CZR CZ half width. 
ay CZSE 4 21-30 kts Hote’ CZ speed error. 
25 CZCE 31-40 deg 0-360 CZ course error. 


4 
26 CZBE 4 41-50 deg 0-360 CZ bearing error. 
4 


27 .GZRE 51-60 mls CZ range error. 


Situation change cards: 


Fields Description ' Fields Description 
1- 5 6-15 

16-20 21-30 

31-35 Input numbers. 36-45 Input values. 

46~50 51-60 

61-65 66-75 


etbid. 
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Table VI. Target Track Summary: Decimal Entries. 


- No. Name Field Units Limits Description 
1 CST(I) 1-10 deg 0-360 Base target course, leg I. 
" 2 PT(I) 11-20 kts pos Base target speed, leg I. 


3 TT(I) 21~30 hrs pos Terminal time, leg I. 


4 =PINI(I) 31-40 0-1 Intelligence detection proba- 
bility, leg I. 

5  SPAX(I) 41-50 mis Surveillance position error 

6 SPAY(I) 51-60 mis parameters, leg I. 


Table VII. Intelligence Data Summary: Decimal Entries. 
No. Name Field Units Limits Description 


1 TI(J) 1-10 hrs pos The time of intelligence 
detection J. 


2 RINP(J) 11-20 mis Intelligence range error, 
detection J. 


3 BINP(J) 21-30 deg 0-360 Intelligence position bear- 
ing from target, detection J. 


4 CINP(J) 31-40 deg 0-360 Intelligence course estimate, 
; detection J. 


5 SINP(J) 41-50 kts pos Intelligence speed estimate, 
detection Jd. 


Table VIII. Terminal Card: Integer Entry. 
Field Description 


ta> Booo9 
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B. DATA DECK ASSEMBLY 

The data deck is composed of a set of environmental 
cards, first run data and the input cards for subsequent 
runs. The environmental set as described in Chapter V, 
Section A-l is placed in front of the deck. The first run 
Gata includes the run identification cards, the option data 
card, the situation data cards and the target track data. 
If intelligence times are to be input, intelligence data 
cards must also be included. The minimum input for subse- 
quent runs is comprised of the run identification cards. 
In addition, there may be an option card or change cards, 
situation cards or change cards, a revised target track, 


and revised intelligence data. Table IX is an outline of 


a typical data set assembly. Figure 19 shows the complete 
data deck assembly. Figure 20 is the example data deck 
used to generate Figures 13 through 16. | 

Job control cards: The exact job control cards are a 
function of the system being used. The job control cards 
required by the IBM system 360 installation at the Naval 
Postgraduate School are listed in Table X. Part A of Table 
X shows the cards required to compile and execute the VIM 
FORTRAN source and data. VIM is currently held in disc 
storage at the NPS computer facility. Part B of Table X 


lists the cards required for use of the stored program. 
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Table Ix. Data Set Assembly 
é Environmental data 
- Title card 
Attacker noise curves 
Target noise curves 
Propagation loss curve 


First run data 
Run identification cards 
Option data card 
Situation data cards 
Target track cards 
Intelligence data cards (optional) 


subsequent runs : 
Run identification cards 
Option data card (optional) 
Option change card (optional) 
Situation data cards (optional) 
Situation change cards (optional) 
Target track cards (optional) 
Intelligence data cards (optional) 


Terminal card 
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Figure 19. Data Deck Assembly for a Series of Two Runs 
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Table X. Job Control Cards 


= A.. Compile and execute 
ie Job card 
te EXEC FORTCLG, REGION.GO=110K 
//FORT.SYSIN DD * 


SOURCE DECK 


/* 
//GO.FTO6FOOl] DD SYSOUT=A,SPACE=(CYL,(1,1)) 
//G0.SYSIN pp # | 
DATA DECK 
/* = 


B. Disc storage 


ta Job card 
//JOBLIB DD DSNAME=S0802 ,KILSUB,DISP=SHR, 
// VOLUME=SER=MARY , UNIT=2314 
//STEPI. EXEC PGM=VIM,REGION=110K 
//FTO6FO01 DD SYSOUT=A ,STACE=(CYL,(1,1)) 
//FT05001 ‘Dp #* 

DATA DECK 
/* 
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VI. VIM EXPERIMENTATION 


The purpose of this chapter is to describe the use of 
VIM in general terms. Section A provides a discussion of 
problem formulation and the relation to the problem in qual- 
itative terms of the input parameters. Section B is an 
illustration of the use of VIM as applied to a typical 


problem. 


A. FORMULATION OF THE EXPERIMENT 

Formulating a simulation experiment involves defining 
an objective, developing a@ scenario compatible with the 
model, enumerating the output data cesired and assigning 
values to the input parameters compatible with the objective, 
the scenario and the desired output. 

Objective: The goals of the experiment must be clearly 
defined in terms of the problem faced and the measure of 
effectiveness to be used. For example, the’ measure of 
effectiveness might be cumulative detection probability or, 
number of counterdetections. 

Scenario: The scenario must relate the weep ien to the 
model. The key elements of the VIM scenario are the trans- 
itting target, a waiting attack submarine, a remote intelli- 
gence facility, the sonar environment, reaction of the 
target to counterdetections, and the fact that the response 


of the attacker to a detection opportunity is to record the 
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probability. of detection and continue its search plan. 

Since the geometry of VIM is similar to that encountered 

in a variety of situations, the elements of the problem 
must be carefully identified with the elements of the model. 
Frequently assumptions must be made regarding the problem 
to satisfy the requirements of the model. Each of these 
assumptions must be carefully justified in light of the 
potential effect on the outcome. 

Output data: Certain data is compiled by VIM automati- 
cally, such as overall detection probability, maximum 
detection aHopabiity. counterdetection data, times associ- 
ated with maximum probability of detection and times associ- 
ated with convergence zone events. Easily obtained from 
the trial summaries is the distribution of detection proba- 
bilities over time for each trial. Other data is available 
from the detailed trial history. 

Input parameters: The variables controlled by input 
were selected to provide a maximum freedom in fitting a 
scenario to VIM. All may have a significant effect on the 
outcome and the values entered heute be the result of 
sivergt edsecuen as well as a clear understanding of the 
input description. 

The environmental data should reflect the best available 
knowledge of the sound propagation profile of the area under 
consideration and the noise characteristics of the submarines 
being studied. Sadaiied of the imprecise nature of the ocean 


environment and the variation between ships of the same 
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class, a parametric analysis of this data group may be in 
order. 

Target leg changes in VIM are intended to interact with 
the intelligence facility as distinguished from interaction 
with the attacker. Changes which would not be discernable 
to the assumed intelligence platform shovld not be included. 
Particularly in the case of intelligence course and speed 
computed from position estimates, frequent leg changes with 
only minor modifications’ will result in unrealistically 
poor intelligence estimates. In general, target track 
should be as simple and direct as the nature of the study 
permits. If random tracks are used it should be only after 
it has been determined that the increased variation is 
necessary to adequately reflect the system being modeled. 

Random intelligence times may be required to simulate 
the rate at which a surveillance facility generates intel- 
ligence estimates on transiting targets. The parameters of 
the random intelligence time option are a time interval 
TINT, SD-13, and the probability PINT, TD-4, that a surveil- 
lance facility will generate an estimate during the interval. 
The intelligence rate is matched by the ratio of PINT:TINT, 
where PINT is between zero and one. If TINT is very small 
a poisson process is approximated. As TINT is increased 
a limit on the maximum frequency of intelligence estimates 
is established with the time between the estimates taking 
on a more regular pattern. For example, assume that it 


takes a submarine six hours to transit an intelligence 
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window, with an average of 3/4 estimates each transit and 

a maximum of 2 estimates on any transit. The pairs (3/4,6), 
(3/8,3) and (3/24,1), where the first number represents 

PINT and the second TINT, would all satisfy the requirement 
of 3/4 estimates per transit. Only the second pair would 
satisfy the limitation that at most 2 estimates could occur 
on any transit. 

Input intelligence times taken from empirical observa- 
tions can be used to validate the model and assist in adjust- 
ing inputs so that exercise results may be extrapolated to 
a more general situation with a degree of confidence. If 
the problem permits this, it is highly desirable. 

The nature of the convergence zone inputs is dictated 
by the problem. If convergence zone intelligence is to be 
generated the error parameters should reflect the errors 
which would be expected in the submarine versus submarine 
tracking problem at convergence zone ranges. 

The number of trials per run is a fuction of the desired 
level of statistical significance. If the parameter to be 
measured is one of the statistics included in the run sum- 
mary, the average over the runs and the sample variance will 
be listed. If the number of trials, N, is large it may be 
assumed that the average value is a normal random variable 
with variance equal to V/N, where V is the listed sample 
variation. Let ¢ be the level of confidence with which the 
results from two runs are to be declared to be from differ- 


ent distributions; i.e., that the two averages represent a 
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different result based on enecaiple pireneheet as distin- 
guished from a random difference. Let Z(c) represent the 
number of standard deviations from the mean of a standard 
normal distribution within which the expected fraction of 
realizations is c. The minimum separation between observed 


averages which may be declared different is 
d = + 2(c)V2V/N 

and if d is specified, the required number of trials is 
N = 2VZ*(c)/d? 


Verification: Before generating the experimental data, 
it is desirable to verify that the program is responding 
properly to the input data. The trial history is an excel- 
lent tool for this task. A sample trial can be plotted 
based on the information provided with each event and the 
pilot examined for signs of program and data errors or 
misinterpretations. 

Validation: Whenever possible the results of the simula- 
tion should be checked with empirical data to determine if 
essential differences exist which might require reformula- 
tion of the problem or if there is sufficient agreement to 
justify confidence in the results where no empirical data 
exists. As pointed out above, the various forms of intelli- 
gence input and computation are particularly well suited 


to this task. 
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B. VIM EXAMPLE 

This typical problem formulation was developed to assist 
the user in becoming familiar with VIM. Initially the 
problem is outlined and the scenario is defined. The control 
parameters are isolated and input data requirements are 
outlined. Finally, the experimental plan is described and 
the results presented. 

Problem: An attack submarine has been assigned the 
task of intercepting at close range a transiting target sub- 
marine. The commander of the attack submarine desires to 
use tactics which will maximize his detection probability. 
The decision variables under the control of the attack sub- 


marine are listed in Table XI. 


Table XI. Example Problem Control Variables 
A list of the variables under the control of the attack 
submarine commander, along with their program names, input 


reference numbers and ranges of permissible values. 


Variable , Name, Ref. Range 

Attacker transit speed SA, SD-9 15-25 knots 
Attacker search speed SS, SD-8 8-16 knots 
Communications interval CINT, SD-2 234 ewes Le neurs 
Distance of attacker YTI, SD-4 0-2000 miles 


from northern boundary 
of corridor 


Duration of attacker Sl, SD=12: ; 0-12 hours 
learch leg 

Speed estimate safety . SAF, SD-19 0-8 knots 
factor 
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seenario: <A potential enemy must use an ocean lane 400 
miles wide and 800 miles long for its submarines to transit 
from home port to their operating stations. During a trans- 
it of the lane a shore based intelligence facility maintains 
intermittent contact with the transitting submarine and 
makes periodic position, course and speed estimates. Upon 
leaving the transit lane the target may head for any one of 
several operating areas. The surveillance system is unable 
to estimate the new course. It is desired that an attack 
submarine be sent to the area to establish sonar contact 
with a transitting target submarine. The political climate 
is such that if the attack submarine is counterdetected the 
target will neither evade nor use weapons. 

Data requirements: The key data in this problem include 
the environmental data, target track data, communications 
delay, the intelligence parameters and number of trials per 
run. The environmental data (including the related situation 
data inputs) is typical of the ships and ocean area being 
considered. No convergence zone is present. 

The target track started at a point north of the attacker 
and two hundred miles to the east and headed due south at 
sixteen knots, the assumed transit speed of the target. 
After 800 miles on the southerly heading, the target turned 
east at the same speed. By this means the 400 mile wide 
corrider was depicted with the attacker on.the centerline of 
the lane and the target proceeding along the boundary. This 


represented the worst case from the point of view of the 
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attacker. The 90 degree turn after 800 miles represented 
a turn to an unknown destination. Surveillance detection 
probability was set to zero on this leg. 

The communications delay was that expected for the sur- 
veillance system. 

Observations from the area of interest indicate that the 
intelligence facility expects four intelligence estimates 
during the course of a target transit with no more than 
three in any 12 hour period. This was approximated by set- 
ting the probability of a surveillance detection during a 
four hour period to one third. It was further assumed that 
the position estimate provided by the attacker would fall 
uniformly within a 40 mile square centered at- the actual 
target position, with course and speed estimates based on 
successive position estimates and speed error limited to 
three knots. 

The number of trials for each run was selected to dis- 
tinguish between average probability estimates differing by 
.03 with 70% confidnnee. A test run of 25 trials was made 
which indicated the variance for papeaiien probability was 


.16. The number of trials required per run was 
N= 2 27(.7)V/a? = 2x(1.05)?x0.16/(0.03)? = 400. 


Since VIM is capable of only one hundred trials per run it 
was necessary to conduct each run four times. 
Experimental plan: A six variable search plan was 


developed to find the combination of variables over the 
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permissible range which would produce the maximum detection 
probability. It was assumed that the probability of detec- 
tion is a smooth function of the six parameters and that a 
single relative maximum exists over the range of the study. 
The variables were searched one at a time for a maximum 
resultant detection probability according to the following 
plan: 

1. A base value and step value were arbitrarily as- 
signed to each variable. 

2. A run (400 trials) was made with each variable at 
its base value. 

3. A run with five variables assigned base values and 
one altered by its step value was conducted for 
each of the variables. 

4, The variable which demonstrated the greatest effect 
on detection probability was identified. A simul- 
taneous search of up to eight values of the variable 
was conducted to estimate the value which maximized 
detection probability. The maximizing value became 
the base for this variable. 

Ds gine the revised base steps three and four were 


repeated until no further improvement was indicated. 


Experimental Results: In the discussion that follows, 
the data from each series of runs is summarized in tables, 
analyzed in view of determining maximum detection probabili- 


ty and the following series specified. A series refers to 
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a set of several runs which are expected to provide the 
answer to the problem of which variables may be adjusted to 
improve detection probability, or what value of a variable 
provides the maximum détection probability for given values 
of the other variables. 

series I was a search from the base case for the vari- 
able which might produce the greatest improvement. Series 
I results are shown in Table XII which lists the base value 
for each variable, the increment by which each was varied, 
the average probability achieved for each variable at its 
incremented value, and the difference of the incremented 
results from the base result. Because of the large number 
of trials for each run (400), the average probability may 
be assumed to be a normally distributed random variable with 
variance V = .16/400 = .0004. The variance of the differ- 
ence of two such random variables is .0008 with standard 
deviation of .0283. The probability that two of the random 
variables with the same mean would differ by .03 is less 
than .3, and the probability of a difference of .06 is less 
than .001. With better than 99% confidence, then, it can 
be inferred that the expected detection probability with 
transit speed set at 22 knots is greater than the expected 
result with transit speed of 20 knots. Series II was a 
search of values of transit speed to estimate the maximiz- 
ing value. 

In Table XIII is listed the results of varying attacker 


transit speed in one knot increments from 20 to 25 knots. 
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Table XII. Series I Results 


Variable CINT =: YTT SS SA SI SAF Base 
Case 

Base Value 4 400 13 20 cast 4 

Increment ~2 —50 +2 | +2 +.23 -.5 

Probability 4575 4425 870.5375. 4750~— 4750 ~—-. 70 

Difference ~.0125 ~.0275 0 +.0675 ~.0325 +.005 


Table XITI. Series II Results 
Search of attacker transit speed. 


SA 20% eal _ 2et 23 wat en 


Probability .470 .460 519 4825 «5250 .5850 


* Series I results. 


Figure 21 is a plot of the results with a least squares 
regression line indicating the trend. Since the results 
seem to be increasing with transit speed and transit speed 
is constrained less than 25 knots, the maximizing value is 
taken to be 25 knots which is used as the new base value. 
Series III was a check of all variables one at a time using 
the revised base. 

Table XIV lists the revised base values and increment 
values for each variable along with the resultant probabili- 
ties and differences. In this case communications interval, 
CINT, showed the greatest deviation, while search speed, SS, 
initial position, YTI, and estimation safety factor, SAF, 


all showed significant deviation. 
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Figure 21. Series II Data with Least Squares Regression 
Line. 
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Table XIV. Series III Results 


Variable CINT SUE. ss SI SAF ALL 
Base Value i 400 13 2d) | 4 BASE 
Increment +4 +100 -3 +.23 -1.5 VALUES 
Probability 04575 25425 4850 5650 5375 .5875% 
Difference -~.1275 -.0425 -.10 -.02 -.0475 


* Series II results. 


Series IV was a search of the permissible communications 
intervals, CINT, for the maximizing value. Table XV shows 
the results of this rearch. The values achieved with CINT 
equal to 6 or 8 are obviously lower than those for CINT 
equal to 2 or 4, but there is insufficient information to 
distinguish between 2 and 4 hour communications intervals, 
and the base value was left at 4 hours. 

Series V commenced a search of attacker search speed 
and interval. The base interval represents a search leg of 
30 miles at 13 knots. It is reasonable to suppose that a 


coordinated adjustment of search interval and search speed 


Table XV. Series IV Results 


Search of communications interval 


CINT 8 6 . 4 2 


Probability 45.75% 50.75 56.37 57.0 


* Series III results. 


@ Average of Series II and Series IV results. 


126 





might have greater effect on detection probability than 
either variable leg “itself. Series V examines search speed 
alone. Series VI examines search interval alone. Series 
VII and VIIT look at the joint variation of search speed 
and search interval. Table XVI lists the results of the 


search speed experiments. Although there is a peak at 13 


Table XVI. Series V Results 


Search over values of search speed 


SS 10 lee 13 14 16 18 20 


Probability 51.25 53. 56.37% SH 52° 46.25 46.5 


® Average of Series II and Series IV results. 


knots, there is little statistical difference over the 
range 12 to 14 knots. Thirteen knots was retained as the 
value of search speed. 

Table XVII lists the results of the search interval 
experiment with search speed fixed at 13 knots. Figure 22 
is a plot of the results. The table and figure indicate 
a degree of stability over the range of 1/2 to 3 hours. 

The average result from all runs for which SI is no greater 
than 3 hours is shown in Figure 23. There is no reason to 
reject the hypothesis that all results within this band 

were from the same distribution with the mean equal to the 
sample average. When the average is taken over all runs for 


which the value of SI is no greater than 5, neither the 
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Table XVII. Series VI Results 
Search over values of Search Interval. 
SI 25 nies 135 area: 2.0 2.25 


Probability 54.75 Dad Dkeeo 5G TD BN .25 53829 


St 2/3) 2.5 2.54 2.75 3.0 3.5 
Probability 56.37* 51.5 56.5 Bde 5? ‘b3e5 49 


SI 4 5 8 


Probability 48.5 49.5 43.0 


* Average of Series II and IV results. 


three largest nor the three smallest values within the band 
could reasonably be assumed to come from a normal distribu- 
tion with the revised average value. The conclusion is that 
the results in this case are stable for search intervals of 
up to three hours, but detection probability falls off 
thereafter, 

Series VII adjusts search speed and search interval 
jointly with search interval decreasing by .23 hours for 
each increase of 1 knot in search speed. Table XVIII lists 
the results of this first joint experiment. As in Series V 
stability is indicated for search speed in the range of le 
to 14 knots with no indication of the effect of search 
interval. 

In Series VIII, SI increases .23 hours with each knot 
increase in search speed. Table XIX lists the results of 


Series VIII. Figure 23 is a graph of Series VIII data with 
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Figure 22. Series VI Data with Partial Sample Averages 


Shown. 
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Table XVIII. Series VII Results 


Joint Search over search speed and search interval 


ss 12 13 14 15 16 20 
SI 2.54 2 31: OOF 1.84 le cae .69 
Probability 5350 .5637% 25375 -4975  .5075 4625 


* Average of Series II and Series IV results. 


Table’ ATX. series VIII Results 


Joint search over search speed and search interval. 


« 


Ss 6 7 8 9 - 10 Ld 
SI Pa 93 Te 6 1.39 1.62 1.85 


Probability 5475 e5675 5850 .5850 45.00 55650 


Ss 12 13 14 5 16 20 
SI 2.07 cares vil 254 ar a 3.0 392 
Probability -5750 5650". .5050 4850 -4975  .3925 


* Average of Series II, IV and VIII 


a least squares quadratic regression curve superimposed on 
the data. Under the assumption that the quadratic curve 
represents the underlying relationship between search speed, 
search interval: and detection smibantdey. the maximizing 
values are 8.35 knots and 1.24 hours with an expected detec- 
tion probability of .5792. These values were entered in 

the base case. | 


Series IX is an experiment to find the variable most 


likely to improve the current base. Table XX lists the base 
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Figure 23. Series VIII Results with Least Squares Quadratic 


Regression Curve. 
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Table XX. Series IX Results 


Variable CINT YTL Ste) SA. si SAF ALL 
Base Value 4 HOO 8,35 25 1.24 4 BASE 
Increment -2 -200 42.35 -2.35 -2 +.66 -.66 -2 VALUES 


Probability .6175 .520  .545 5551/5 eS. “b20 25525 .5825 .5792% 


Difference +.0383 -.0592 -.0342 -.0617 .0317 -.0592 -.0267 +.0033 .5792% 
* Least squares estimate. 


, 


values, increments and results for this series. As in 
Series III, the communications interval offers a significant 
improvement. Since the value of communications interval for 
Series IX was its lower bound, no further searching was 
required, and the base was revised to include CINT=2. 

Series X was a search over all variables to find direc-~ 
tions of possible improvement and if no improvement is 
indicated, the range of stability of the solution. Table 
XXI shows the results from Series X. No improvements were 
indicated from variation of any of the parameters. The 
final solution is summarized in Table XXII. The solution 
appears very insensitive to variation of distance from the 
northern boundary, YTI, over the range of 200 to 1000 miles. 
Although most results in this range were greater than the 
base result, the differences were small and did not indicate 
a pattern. The loss resulting from increasing communications 
interval, CINT, from 2 to 4 hours is not statistically signi- 


ficant, while the next step to 6 hours causes a notable 
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Table XXI. Series X Results 


. "Variable YTI YTI YTI YTI YEI YTI 
Value 200 .300 500 600 700 800 
Probability  .580 .590 .6050 45725 .6050 5825 
Difference .0019 = =.0119 0269 ~.0056 .0269 0044 
Variable gle YUL SAF SAF SAF SAF 
Value 1000 1500 i 2 3 a 
Probability 5875 5150 5475 540 «550 25125 
Difference .0094 —-.0631 -.0306 -.0381 -.0181 -.0656 
Variable St/sS S81/SS . SA CINT CINT BASE 
Value 1.16/8 1.32/7.7 24 y 6 CASE 
Probability 25125 25325 ~5375 565 ,5 10 .5781* 
Difference ~.0656 -.0456 ~+.0406 ~-.0131 ~.0131 


* The average value over 1200 trials. 


Table XXII. Final Results 


Variable Solution Range of 
Value Insensitivity 
CINT 2 oayd 
pga 400 200-1000 | 
SS 8.35 , 
SA/SI 25/1.24 % 
SAF ui 3-4 


* The results seem highly sensitive to these parameters. 
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reduction in detection probability. The speed estimate 
safety factor, SAF, below 3 knots or greater than 4 knots 
is associated with a decrease in the final estimate. Mov- 
ing attacker transit speed, SA, away from its upper bound 
or adjusting search speed, SS, and search interval, SI, in 
either direction indicate lower detection probabilities. 
The validity of the results should be tested against 
some of the data estimates. Of particular interest are the 
surveillance parameters regarding the accuracy of the in- 
telligence, the frequency of the surveillance detections 
and communications delay, the speed of the target, and the 
environmental data. It is likely that under differing 
circumstances, the tactics should be altered from those 


derived here. 
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APPENDIX A. LAYER DEPTH ANALYSIS 


This appendix deals with the inclusion of sound velocity 
profile (i.e., layer depth) information in the analysis of 
a problem. The simple case in which only two operating 
depths are available is treated first. The method presented 
is then generalized to the multiple layer depth problen. 

VIM was not designed’ to permit the submarines to change 
depth to take advantage of various sound conditions. Since 
a decision to change depth is a rational choice based on 
the tactical situation, treating it as a random variable 
would complicate an experiment without increasing the relia- 
bility of the results. Further, VIM is not sufficiently 
sophisticated to make depth change decisions on a rational 
basis. 

Layer depth data includes a propagation loss curve 
associated with each of the operating depth combinations. 
The objective of this analysis is to synthesize a single 
propagation loss table which reflects the expected detection 
ranges resulting from optimal operating depth policies by 
both submarines. By considering the depth selection problem 
as a two person, zero sum game [3] this objective can be 
realized. 

In the two person, zero sum game format the attacker 
and target are both faced with deciding the proportion of 


time to spend at each operating depth. A payoff, expected 
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detection range, is associated with cach depth depending on 
the operating depth of the other submarine. Finally, a 
desirable payoff from the viewpoint of the attacker is a 
long detection range, whereas this payoff is undesirable 
for the target. 

1. One Layer Problem 

Assume that at some given depth there is a layer 
which either ship can choose to be above or below. From 
the target radiated noise curve and the attacker self noise 
curve determine the noise levels for representative speeds 
and compute the figure of merit for detection of the target 
by the attacker. Using the figure of merit and the appro- 
priate propagation loss curves determine the maximum ranges 
of zero signal excess for each combination of operating 
depths. The problem can be summarized in matrix form as 
shown in Table A-l. 

The value of the game matrix is the expected payoff 
when both submarines are operating to their best advantage. 
This is the detection range which will be matched with the 
figure of merit to form the input propagation loss tables. 

The first step in finding the value of the game 
matrix of Table A-1 is to test for dominance of one operat- 
ing depth over another for either ship. If the target is 
always better off operating at depth 1 than if it used 
depth 2, regardless of the depth of the attacker, then tar- 
get depth 1 dominates target depth 2 and the column asso- 


clated with depth 2 may be deleted from the analysis. 
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Table A-1 


ATTACKER 


Table A-2 


One layer, two operating depth game matrix. 
Entries represent expected detection ranges 
when the submarines are at the depths in 


indicated. 
TARGET 

Operating Operating 
depth (1) depth (2) 
Above the Below the 
layer layer 

Operating 

depth (1) 

Above the 

layer 

Operating 

depth (2) 

Below the 

layer 





Example of two layer, three operating depth 

game matrix. Entries are the expected detection 
ranges when the submarines are at the depths 
indicated. 


TARGET 


ATTACKER 
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In Table A-1, R,, would be less ae R,, and Hay would be 
less than R,. Having eliminated one possible decision for 
the target, the attacker may now choose the operating depth 
corresponding to the larger of Ry, and Ry which is the 
value of the game. 

For example, if the actual matrix of detection 
ranges were 

Target 


Depth 1 2 


. FE 
Attacker 
2 pepe 


then the target would always select depth 1 and the attacker 
would remain in depth 2. Expected detection range would be 
le miles. 


Similarly, if the matrix were 


Target 
Depth 1 2 


1 8 10 


Attacker 


then Seeaenens depth 2 is always better than depth 1, and 
the target would choose depth 2 resulting in an expected 
detection range of 11 miles. 

If it is established that no dominance exists between 
the rows or columns of the matrix, the solution for the 
expected detection range, R, is found in reference [4] and 


is given by 
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Ru [Ra - Re | + Ra [Ru - Riz | 


[Ra = a | + [Rn - Riz | 


as Multiple Layer Problem 

‘Where more than two operating depths are permitted 
the nature of the problem is unchanged but the mathematics 
required for the solution are more complex. The detection 
matrix for each combination of depths is now N x N, where 
there are N possible operating depths. For the purposes 
of this discussion let N=3 and assume Table A-2 is the ex- 
pected detection range matrix for a given figure of merit. 
As in the one layer, two operating depth case the first 
step is to test for dominance. Check to see if each of the 
entries from a row or column are preferred to the corres- 
ponding values of another row or column. If so, delete 
the dominated row or column from the analysis. In this way 
it may be possible to reduce the problem to simpler form. 
In the example, Table A-2, there is no dominance. 

After removing dominated rows or columns, test the 
matrix for a saddle point. For each protagonist list the 
worst possible outcome from each possible decision. By 
selecting the operating depth which results in the least 
undesirable of these outcomes each submarine is guaranteed 
that the expected detection range will never be worse than 
that value regardless of the actions of the other. If the 
range so chosen is the same for both submarines their deci-~ 
Sions will actually result in that expected value and a 


saddle point has been found. 
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In the example, Table A-2, let A, = min Rey and 


T., = max Req: 
J a J 


A; = min (10, 12, 14) = 10 
Ao Som Cb ASB) 
A; = min (13, 17, 12) = 12 
T,; = max (10; 25, 13) 5 15 


T. = max (12, 13, 17) = 14 








T, = max (14, 9, 12) = 14 


Let A*® 


> 
HT 


max max (10, 9, 12) = 12 


ae 
min (15, 17, 14) = 14. 


and T*® = min T 
J 
Since A®* # T#, the example does not have a saddle point. 


j 


J 


Whether or not the above operations produced any 
simplifications, the problem can be formulated as a linear 


program. The general form is 


maximize R 


subject to R - Ry Xacee  oR yy ky €< 0, k=l, aN 
xit 1 tXy = 1 
Ry Xpoe reo Xy 2 0 


where Xy represents the proportion of time the attacker 
spends at depth i. The actual values of x are not impor- 
tant for this analysis. The value of R will be the expected 
detection range at the given figure of merit. In the exam- 


ple, Table A-2, the problem can be written 


140 





maximize R 


subject to R ~- 10x, — 15x, - 13x, <.0 


fo Dis os ee ee RS SO 
R - Lux, - 9x, - 12x, < 0 
xX, + Kot Re 


WwW 
© 


Ry X)5 Xoo %3 


Solution to this problem by the simplex method is 
explained in reference [4] and programs to carry out the 
computations are available at many computer facilities. 

By selecting seven values of figure of merit near 
those values expected during the execution of the experiment, 
a propagation loss table will be developed which reflects 


optimal submarine operating depth selection. 
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APPENDIX B. PROGRAM LISTING 


1. Alphabetical subroutine index 


Program name 


MOLD yao Sapa eS a eee ee eer et 17e 
CPAT ----------- ne ee ene nee eee 168 
COURSE --------------- 3 -- n-ne e- 184 
CURVE --------- mie rn ne en nn en Hee 181 
inp e oe ea Geese tte at etal w-o----- 170 
CZINT ----------~---~-----~---~----~--- 178 
DEG werner ee en ee er re eee eee 184 
DETECT -----~--~-~-----~-------------- +--+ 169 
ELIST -----~------~----++~~-- ene 178 
PIND -e3 nen n ene nner ren nnn nnn en 167 
INOUT -----~---------~----~~-~------~---~-------- 153 
INOQUT2 --------~----- He --- 154 
INOQUT3 ----~------------------ 156 
INOUTY ------------------------------~----- 157 
INOQUT5 ------~----------------------------- 160 
INOUT6 --~--------------------------~+------- 161 
TNOUTY oneesene eee eee ee es 161 
INOUT8 ----~----~---------~------------------- 157 
INTDET --------~-----~----~------------~------- 175 
; PUG: eevac nee en ieee Joel ees et ee 174 
MAIN -sen-e ene nn nase nnn nnn nn enn nnn nem 144 
NORMAL ---------------~------------------------ 186 





Alphabetical index to Appendix B (cont.) 


Program name 
ORDER ~--------3-- +e ee ne en eee 182 
. ORDER seecieeeo eee cee eee 182 
ORDER3 --~-~--~-~----~- + -~-- -- 182 
PDET ------ n-ne ee en een eee 172 
PIM ~-~-~------~------~-------~~----------------- 183 
PROINT --~---~---------~-—-----~---~---~----~---- 176 
RAD ~~--~-~--~~~---~ + ne ne eee 184 
RAND ~-~---~--~---~-~--~-~~---~--~- ~~~ -- + 186 
REVENT ~~--~—-~-~-.--~------~-- +--+ - ~~ 3 185 
BiG eg ssa ae ee Soe a ects cee 185 
FDL Be ice a he dal ok See chs aris 180 
0S gene ere na re GER ROE HEE rears geOgE TEP SRROTER ef te ep ee 187 
SU weet Seco ee ei ee So ee 180 
UO tee te ee ie ete 167 
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